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INTRODUCTION 


The  Center  for  Thin  Film  Studies  was  formed  in  1986  under  grant  from  the 
Air  Force  Office  of  Scientific  Research  as  part  of  the  University  Research  Initiative 
Program  (URIP), 

This  report  covers  the  first  year  of  operation  of  the  URI  Thin  Film  Center 
(TFQ,  and  describes  a  diverse  array  of  studies  on  thin-film  materials,  substrates,  and 
their  processing  and  analysis.  Individual  efforts  are  highlighted  in  sections  on 
nucleation  studies,  ion-assisted  deposition,  Rutherford  backscattering  spectrometry. 
Brillouin  scattering,  a  continuum  theory  of  the  evolution  of  structure  in  thin  films,  a 
study  of  polishing  parameters  relevant  to  the  preparation  of  substrates,  and  the  setup 
of  a  characterization  facility  for  the  Center. 


The  work  described  herein  was  performed  under  Contract  ti  F49620-86-C- 
0123P0000J 
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TASK  t 


I 

I 

THIN-FILM  NUCLEATION  STUDIES: 

CONVENTIONAL  AND  ATOMIC  LAYER  EVAPORATION 

Principal  Investigator:  U.  J.  Gibson 

Project  Goal 

i 

Study  the  growth  of  thin  films  produced  by  different  deposition  methods,  and  | 

correlate  nucleation  behavior  with  the  development  of  film  microstructure.  ! 

I 

Summary  j 

r  ^ 

During  the  first  year  of  the  project  we  investigated  ZnS  made  under  varying 
conditions,  made  preliminary  measurements  of  substrate-nucleation  effects  in  ZnS 
deposited  on  Si,  and  designed  and  assembled  (minus  one  piece)  a  preparation  chamber 
to  allow  detailed  nucleation  studies  under  ultra-high  vacuum  (UHV)  conditions.  The 
experiments  performed,  equipment  status,  and  plans  for  the  coming  year  are  discussed 
below. 

Work  Accomplished 

Film  Growth 

Thin  films  of  ZnS  have  been  prepared  by  resistive  evaporation  in  a 
conventional  system  under  varying  deposition  conditions,  and  by  molecular  beam 
evaporation  (MBE)  onto  both  single  crystal  and  amorphous  substrates.  These  films 
have  been  studied  with  a  wide  variety  of  techniques,  including  x-ray  diffraction, 
waveguide  loss  measurements,  second-harmonic  generation,  and  Brillouin  scattering. 

In  a  diffusion-pumped  system,  we  found  wide  variations  in  the  structure  of  ZnS 
deposited  onto  silica  substrates  under  different  conditions.  Some  of  these  changes 
(e.g.,  increased  grain  size  with  increased  substrate  temperature)  can  be  readily 
predicted  on  the  basis  of  simple  film-growth  models.  In  addition,  we  found  strong 


variations  in  optical  properties  and  crystal  structure  of  the  films  with  changes  of  the 
residual  gas  composition.  In  particular,  we  found  that  residual  oxygen  in  the  system 
could  reduce  absorptive  losses  of  the  material,  while  changes  in  the  H20  partial 
pressure  produced  dramatic  changes  in  the  preferred  orientation  and  size  of 
crystallites  in  films  deposited  onto  ambient  temperature  substrates.  Changing  the  H20 
partial  pressure  from  3x10^  to  4xlO~T  torr  changed  the  preferred  orientation  from 
(111)  to  (220)  parallel  to  the  substrate,  and  the  sizes  of  the  dominant-orientation  grains 
from  about  130  A  to  greater  than  250  A.  These  changes  strongly  affect  the 
waveguide  losses  of  these  materials,1  as  shown  in  Fig.  1. 


Figure  1.  Propagation  loss  dependence  on  Oj  partial  pressure  during  deposition.  Squares 
and  diamonds  represent  H20  partial  pressures  of  4x10' 7  torr  and  3x10'*  torr  respectively. 

MBE  growth  of  ZnS  onto  Si  substrates  resulted  in  high-quality  crystalline  films, 
as  determined  by  x-ray  diffraction.  Recent  Brillouin  measurements  indicate  the 
presence  of  an  acoustical  mode  analogous  to  that  found  in  ZnSe  (refer  to  the 
discussion  in  the  Task  3  report). 


Second  harmonic  (SH)  measurements  on  ZnS  films  are  also  underway  to  assess 
the  effect  of  vapor  impingement  angle  on  the  structure  of  the  films.  Previous  resuits 
on  alumina  films  indicate  that  SH  is  an  extremely  sensitive  function  of  film 
anisotropy.1 


Nucleation  Studies 


One  of  the  goals  of  this  project  is  investigation  of  atomic  layer  evaporation 
(ALE),  where  reactive  pulses  are  brought  to  a  substrate  sequentially.  In  the  growth 
of  ZnS.  possible  reactants  are  Zn  and  thermally  decomposed  HjS.  Preliminary 
experiments  have  been  performed  to  characterize  the  interaction  of  thin  layers  of  Zn 
deposited  onto  Si  surfaces.  Equivalent  thickness  layers  (as  measured  by  a  quartz 
crystal  monitor)  deposited  onto  different  substrates  gave  widely  varying  x-ray 
photoelectron  spectroscopy  (XPS)  signals,  indicating  different  nucleation  behaviors. 

MBE  layers  of  ZnS  grown  on  etched  and  oxygen-passivated  Si  wafers  indicated 
that  heating  the  wafers  in  the  presence  of  an  HjS  flux  prior  to  deposition  improved 
film  properties,  probably  due  to  a  Si-S  layer  created  via  an  O-S  surface  exchange 
reaction.  Two  sets  of  experiments  were  performed  on  a  thin  Zn  layer  deposited  onto 
Si  prepared  in  this  fashion.  The  first  was  on  the  evolution  of  the  Zn  signal  with 
increasing  temperature,  and  the  second  was  angle-resolved  photoemission  to  determine 
the  nucleation  behavior  of  the  Zn. 

The  first  set  of  experiments  was  used  to  study  the  appropriate  temperatures  for 
ALE  using  elemental  Zn.  Increases  in  the  evaporation  rate  are  expected  at:  a)  the 
Zn-Zn  bond  energy  equivalent,  where  Zn  will  evaporate  from  the  exposed  surface, 
and  b)  the  temperature  equivalent  of  the  Zn-substrate  bond  energy.  We  observed 
slope  variations  at  approximately  275°C  and  550°C,  corresponding  to  these  thresholds. 

Angle-resolved  measurements  were  used  to  test  the  hypothesis  that  sulfur 
treatment  of  the  Si  surface  led  to  formation  of  a  complete  monolayer  as  opposed  to 
clusters  of  Zn.  Although  the  analysis  is  complicated  by  spurious  oxygen  arising  from 
the  transfer  of  the  sample  from  the  deposition  system  to  the  XPS  system,  the  results 
(Fig.  2)  are  convincing  indications  that  the  Zn  does  form  an  unbroken  layer,  suitable 
for  ALE  growth. 

Hardware 

Design  and  acquisition  of  almost  all  the  pieces  of  the  treatment  chamber  to  be 
attached  to  the  XPS  system  have  been  completed,  although  we  still  need  the  transfer 
rod  (ordered  in  February.  1986).  The  configuration  is  shown  in  Fig  3.  An 
ultraviolet  source  is  presently  being  added  to  the  surface  analysis  system. 


Zn/Si  vs.  theta,  assuming  no  oxide  layer. 


□  Experimental  data 

- 100%  coverage 

- 90%  coverage 
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Figure  2.  Comparison  of  theoretical  fit  with  actual  data  for  Zn/Si  (99eV)  area  ratio. 


Collaboration  with  Dr.  Dereniak’s  group  will  allow  us  to  assess  the  possibility  of 
adding  a  scanning  tunnelling  microscope  to  the  system;  a  flange  for  that  purpose  is 
on  the  chamber. 


SAMPLE  xPS  UPS 

PREPARATION  CHAMBER 

CHAMBER  CHAMBER 


URFACE  ANALYSIS  SYSTEM 


Figure  3.  Configuration  of  the  XPS  system  treatment  chamber. 
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Plans  for  the  Coming  Year 

The  plans  for  the  coming  year  center  around  the  new  deposition  chamber, 
which  will  allow  us  to  perform  UHV  experiments  on  the  nucleation  behavior  of  MBE 
and  ALE  films.  Sequential  deposition  and  analysis  steps  will  be  possible  between  the 
two  halves  of  the  system.  Specific  studies  to  be  undertaken  include  ion  activation  of 
the  Si  surface,  and  more  detailed  studies  of  the  HjS-prepared  Si  substrates.  Further 
studies  of  thick  samples  grown  by  MBE  in  the  other  chamber  are  also  planned,  both 
in  collaboration  with  Dr.  Stegeman's  lab,  and  using  the  expanding  capabilities  of  the 
characterization  facility. 
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TASK  2A 


STRUCTURE  MODIFICATION  BY  ION-ASSISTED  DEPOSITION 
Principal  Investigators:  H.  A.  Macleod  and  M.  R.  Jacobson 
Project  Goals 

Investigate  the  use  of  ion-assisted  deposition  (IAD)  for  modification  of  film 
structure,  and  determine  its  utility  for  improving  the  performance  of  optical  films. 

Summary 

Progress  on  the  URIP  contract  has  been  concentrated  in  four  areas,  as  listed 
above  in  the  section  on  project  goals.  The  following  sections  deal  with  each  of  these 
areas  separately,  but  it  should  be  noted  that  there  was  significant  interaction  between 
these  projects;  for  example,  the  IAD  work  relied  heavily  on  RBS  analysis  at  the  Ion 
Beam  Analysis  Facility  and  on  other  characterization  methods  provided  by  the 
Measurement  Laboratory.  We  point  out  here  that  unless  otherwise  noted,  the  coatings 
studied  were  deposited  in  the  Balzers  760  evaporation  system.  This  cryo-pumped  box 
coaler  achieves  base  pressures  of  about  2x10“*  mbar,  and  is  equipped  with  a 
Commonwealth  Scientific  3-cm  Kaufman  hot-cathode  ion  source. 

Work  Accomplished 

Ion- Assisted  Deposition  of  Oxides,  Nitrides,  Oxynitrides.  Fluorides,  and  Metals 

Extensive  IAD  work  has  been  carried  out  on  metal  oxide,  oxynitrides,  nitrides, 
and  fluorides,  and  on  metals,  with  considerable  success.  Parameters  that  have  been 
varied  include  beam  current  and  voltage,  ion  beam  duty  cycle  (continuous  and 
intermittent  bombardment),  and  bombardment  gas  species  (helium,  neon,  argon,  and 
krypton,  as  well  as  oxygen  and  nitrogen).  Summaries  of  the  results  for  these 
materials  follow. 
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Metal  Oxides  -  Alumina  and  Titania 

Bombardment  of  alumina  by  both  oxygen  and  argon  ions  was  conducted  largely 
to  investigate  the  unusual  behavior  noted  in  a  previous  study  by  Ebert.1  Alumina 
exhibits  a  very  high  packing  density  for  oxide  films,  even  when  deposited  on  ambient 
temperature  substrates.  This  is  shown  by  the  lack  of  a  measurable  vacuum-to-air 
shift  in  optical  properties  (comparison  between  properties  immediately  after  deposition 
while  still  under  vacuum  and  after  venting  to  atmospheric  pressure).  The  fact  that 
no  shift  occurred,  even  for  low  index  samples,  indicates  that  alumina  is  dense  enough 
to  preclude  moisture  penetration.  Our  results  also  emphasize  the  importance  of 
oxygen  partial  pressure  on  both  unbombarded  and  bombarded  coatings.  We  noted  a 
strong  increase  in  index  with  LAD,  up  to  1.70  at  3S0  nm,  without  any  correlation  to 
or  change  in  density.  The  films  were  clearly  amorphous  according  to  x-ray 
diffraction  (XRD)  and  transmission  electron  microscope  (TEM)  measurements;  changes 
in  index  are  presumably  caused  by  alterations  in  the  amorphous  network,  rather  than 
in  crystallinity.  Ar  bombardment  increased,  while  O  bombardment  decreased 
absorption;  we  believe  that  O  atoms  can  fill  vacancies  without  adding  impurities. 
Finally,  the  Kaufman  guns  cannot  provide  long-term  O  bombardment  due  to  oxidation 
of  their  filaments.  Cold  cathode  or  filament-free  Kaufman  guns  may  be  preferable. 

Intermittent  oxygen  ion  bombardment  of  titania  films  was  attempted,  to 
determine  whether  it  could  approach  the  increased  refractive  index  and  decreased 
absorption  achieved  in  continuously  bombarded  films.  The  experiment  was 
performed  in  the  Balzers  760  chamber,  with  the  substrates  rotating  through  a  300  V 
ion  beam  from  the  ion  gun  while  accumulating  material  from  thermally  evaporated 
Ti,Oj  tablets.  The  results  were  positive,  with  the  optical  constants  reaching  values 
reported  by  others  for  the  LAD  of  titania. 

Metal  Nitrides  -  Aluminum  Nitride 

We  compared  two  possible  means  of  producing  AIN:  reactive  evaporation  of 
A1  in  an  Nj  atmosphere,  and  IAD  bombardment  of  thermally  evaporated  A1  by  a 
neutralized  N  ion  beam.  The  former  did  not  succeed  in  incorporating  N  into  the 
film,  but  under  the  right  conditions,  the  latter  yielded  stoichiometric  AIN,  or  even 
superstoichiometric  A1NX,  with  x  as  high  as  1.51  We  believe  that  the  N  atoms  were 
implanted,  the  bombardment  stimulated  N2  sorption,  or  that  the  N  atoms  in  the 
plasma  were  somehow  activated  with  respect  to  neutral  N2. 

Films  with  x  <£  1  were  dielectric,  but  had  considerable  visible  absorption. 


uw mp&jV'JL\wJiKmJii\rvvv'j*ur\rw ij\ .-v > -  v -v*  -  *  - j>  - j»  -jt  v  v  v;, 


Annealing  at  500°C  reduced  the  absorption  and  moved  it  toward  the  UV,  presumably 
by  permitting  interstitial  N  to  diffuse  from  the  film  or  O  to  fill  N  vacancies.  A1NX 
films  with  x  as  low  as  0.5  were  also  produced,  but  these  were  two-phase  films,  with 
transparent  and  opaque  crystallites  alternating  throughout  the  film  on  a  scale  of  a 
few  /a n.  For  this  work,  the  3 -cm  gun  was  marginally  able  to  bombard  the  sample; 
the  12-cm  gun  purchased  this  year  with  URIP  funds  will  enable  us  to  further 
explore  this  method. 

Metal  Oxynitrides  -  Aluminum  Oxynitride 


We  explored  the  possibilities  of  tailoring  the  index  of  AIN  by  bombarding 
thermally  evaporated  A1  film  with  N  in  an  02  atmosphere.  This  highly  hybrid 
reactive  IAD  approach  yielded  films  ranging  in  index  from  1.65  to  1.93,  depending 
on  the  ion-gun  current  density.  In  addition,  absorption  was  systematically  lower. 
We  believe  that  oxygen  fills  the  vacancies  left  by  N  bombardment. 

Metal  Fluorides  -  Cryolite  and  Magnesium,  Aluminum,  Lanthanum,  Cerium, 
and  Neodymium  Fluorides 


Despite  the  differences  among  the  many  metal  fluorides  deposited,  our  work 
produced  these  general  conclusions.  First,  IAD  tends  to  sputter  fluorine  atoms 
preferentially,  resulting  in  substoichiometric  films.  Second,  the  resulting  vacancies 
are  often  filled  with  O  or  OH",  which  are  available  from  residual  water  vapor  or 
from  implantation  if  O  ion  bombardment  is  employed.  Visible  absorption  will  be 
minimal  if  enough  O  is  available  to  fill  these  vacancies,  but  some  UV  absorption 
characteristics  of  these  metal  oxides  will  appear.  However,  if  preferential  sputtering 
creates  too  many  vacancies  to  fill,  we  observe  visible  and  UV  absorption. 

Refractive  index  and  packing  density  tend  to  increase  under  IAD.  Other 
related  properties  change  as  well:  vacuum-to-air  shifts,  water  absorption  bands,  and 
tensile  stress  decrease.  While  these  quantities  can  approach  bulk  values  during  more 
intense  IAD,  there  is  a  point  at  which  growing  film  absorption  begins  to  undermine 
these  improvements.  We  have  also  noted  that  crystalline  fluorides  (LaF3,  CeF3,  and 
NdFj)  will  transform  from  fluoride  to  oxyfluoride  phases  when  O  content  reaches  a 
high  level. 

Tungsten  filament  erosion  in  the  ion  gun  measurably  contaminated  the  films, 
according  to  RBS  measurements.  The  contamination  did  not  seem  to  affect  the 
optical  properties. 
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Metals 


The  IAD  of  metals  has  been  studied  much  less  than  the  IAD  of  dielectrics.  In 
the  interests  of  rectifying  this  neglect,  Chang-Kwon  Hwangbo  began  an  investigation 
of  the  IAD  of  Ag.  In  general.  Ag  films  are  deposited  with  significant  tensile  stress. 
IAD  with  Ar  tends  to  reduce  this  stress,  and  even  leave  the  stress  compressive.  The 
effect  appears  to  be  proportional  to  ion  beam  momentum,  rather  than  beam  energy  or 
current  density.  Momentum  also  correlates  with  changes  in  optical  constants,  lattice 
spacing,  electrical  resistivity,  and  incorporation  of  Ar  ions. 

RBS  Analysis  Techniques 

During  the  past  year,  we  have  worked  extensively  with  Prof.  J.  Leavitt, 
director  of  the  Ion  Beam  Analysis  Facility  (IBAF)  at  the  University  of  Arizona 
Physics  Department.  There,  a  large  Van  de  Graaff  accelerator,  described  elsewhere 
in  this  report,  provides  a  well-characterized,  high-current  ion  beam  for  analysis  of 
surfaces.  Our  main  effort  has  been  the  derivation  of  a  correction  for  heavier 
elements  (high-Z  correction)  to  their  absolute  Rutherford  scattering  cross-sections  {a). 
The  conventional  calculation  of  or  ignores  the  effect  of  the  electrons,  simplifying  the 
encounter  to  a  simple  two-body  collision.  This  approximation  holds  very  well  at 
high  collisional  energies  (>  3  MeV)  and/or  low  atomic  number  for  the  bombarded 
atom.  However,  for  high-Z  elements  and  energies  below  3  MeV.  a  common  situation 
in  RBS,  a  must  be  corrected  to  account  for  electron  shielding  of  the  nucleus,  which 
effectively  reduces  the  nuclear  charge  experienced  by  the  incoming  particle. 

The  high-Z  correction  is  not  critical  for  collisions  between  particles  of  similar 
atomic  numbers,  but  becomes  important  when  the  particles  are  separated  by  more 
than  a  row  in  the  periodic  table.  If  the  correction  is  ignored,  the  error  may  exceed 
the  normally  quoted  precision  of  the  technique.  Moreover,  computations  of  atomic 
densities  from  RBS  data  must  include  the  correction  regardless  of  the  magnitude  of 
the  atomic  mass  variations.  We  have  verified  the  importance  of  the  high-Z 
correction  in  both  regards. 

In  a  related  area,  we  produced  an  analysis  program  for  IBM-PC  standard 
microcomputers  to  expedite  RBS  analysis. 


Simulation  of  Thermal  and  Growth  Processes 


Although  explicit  funding  for  this  area  begins  in  October,  1987,  we  have 


carried  out  some  work  with  the  assistance  of  Robert  Sargent,  who  is  supported  by  a 
Kodak  Fellowship.  He  vastly  improved  the  model  code,  and  has  modified  it  to  run 
on  the  MicroVax  computer  here  at  the  Center.  The  results  have  been  applied  to 
models  of  defect  propagation  and  multilayer  deposition.  In  addition,  we  also  began 
theoretical  work  in  thermal  modeling  while  partially  supported  by  this  contract,  and 
organized  an  SPIE  Symposium  on  the  topic  of  computer  simulation  of  the  growth  of 
thin  films  for  the  1987  Annual  Meeting  in  August. 

Characterization  Facility 

In  the  interests  of  assembing  a  central  characterization  facility,  we  have 
successfully  solicited  the  donation  of  a  Cary  2415  spectrophotometer  from  the  Varian 
Corporation.  This  instrument,  along  with  the  upgraded  Fourier  transform  infrared 
(FUR)  spectrophotometer,  will  complement  the  electron  microscopes  and  WYKO 
surface  profilomster  obtained  for  the  characterization  facility. 

Plans  for  the  Coming  Year 

Ion- Assisted  Deposition  of  Oxides,  Nitrides,  Oxynitrides,  Fluorides,  and  Metals 

LAD  will  continue  to  be  the  area  of  concentration  in  the  second  year.  We  plan 
to  work  on  a  number  of  problems  which  were  briefly  investigated  during  the  first 
year: 

*  Commissioning  of  the  12-cm  ion  gun  -  because  of  long  delays  related  to 

bidding  and  Arizona  law,  we  did  not  receive  the  12-cm  gun  and  its  power  supply 

until  September.  We  are  currently  installing  it  in  the  Balzers  760.  and  it  will  be 
available  during  the  coming  months  to  enhance  the  size  and  uniformity  of  our  ion 
beam  for  several  experiments. 

*  EXial  ion  beam  experiments  -  the  arrival  of  the  new  gun  allows  us  to  begin 

experiments  where  the  smaller  gun  can  ion  sputter  source  or  target  material  and  the 
larger  gun  can  bombard  the  growing  film.  Such  dual  km  beam  capability  will 

multiply  the  parameters  under  our  control  during  deposition. 

*  Intermittent  bombardment  -  preliminary  work  indicates  that  intermittent 

bombardment  of  titan ia  yielded  improvements  similar  to  those  of  continuous 

bombardment. 


We  intend  to  extend  this  work  to  other  materials: 

*  Angular  dependence  -  we  will  investigate  whether  the  angle  at  which 
bombarding  particles  arrive  at  the  substrate  affects  the  ultimate  properties  of  the 
coating. 

*  Exotic  gas  bombardment  -  while  argon  is  the  common  gas  employed  for  LAD. 
due  to  its  moderate  mass  and  low  price,  it  is  possible  to  bombard  the  film  with  other 
noble  gases  and  nitrogen  to  permit  variation  of  momentum  independently  of  energy. 
Preliminary  results  with  He,  Ne,  and  Kr  indicate  clear  index  dependences  on 
momentum  rather  than  energy. 

*  Materials:  similar  oxides,  nitrides  and  fluorides  will  be  deposited  this  year. 
For  metal  films,  we  will  embark  on  a  program  of  A1  deposition.  It  is  known  that 
Al.  in  contrast  to  Ag.  is  ion-beam  sputtered  in  a  compressive  state.  As  we  have 
done  for  Ag.  we  will  compute  optical  constants  with  the  surface  plasmon  resonance 
technique,  determine  microstructural  properties  with  XRD.  measure  electrical 
resistivity  with  the  four-point  probe,  and  derive  stoichiometry  with  the  RBS 
technique.  It  has  rlso  been  discovered  recently  that  silver  can  be  efficienlty  oxidized 
in  a  plasma,  so  that  oxygen  bombardment  of  silver  may  be  expected  to  yield  some 
interesting  materials. 

RBS  Analysis  Techniques 

We  will  continue  to  subject  our  samples  to  RBS  analysis  and  to  improve  the 
technique: 

*  Better  elemental  resolution  -  the  ongoing  work  at  the  RBS  facility  to  double 
the  beam  energy  to  over  10  MeV  by  doubly  ionizing  the  projectile  nuclei  should  be 
completed  this  year.  Under  certain  conditions,  this  can  be  used  to  simplify  elemental 
identification  by  reducing  overlap  between  peaks,  producing  higher  resolution  between 
elements  in  thicker  films.  This  improvement  will  be  applied  immediately  to  the 
analysis  program. 

*  Nuclear  reaction  analysis  for  H  profiling  -  we  will  also  participate  in 
employing  the  nuclear  reaction  analysis  (NRA)  capability  that  the  IBAF  intends  to 
complete  this  year.  To  probe  the  distribution  of  hydrogen  with  depth,  the  most 
difficult  atom  to  observe  for  most  techniques,  a  1SN  projectile  nucleus  is  accelerated 
to  an  energy  of  6.4  MeV  to  10  MeV.  The  l5N  nuclei  react  with  hydrogen  atoms  and 
produce  a  '*C  nucleus,  an  a  particle,  and  a  4.43  MeV  7-ray.  The  resonance  can  be 


easily  observed.  The  method  is  especially  valuable  in  profiling  water  in  surface 
layers.  Since  H*  molecules  are  rare  in  the  atmosphere,  we  can  attribute  observed  H 
either  to  H,0  or  to  OH~  which  has  been  physisorbed  or  chemisorbed  to  the  film 
surface.  By  comparing  H  and  O  results,  we  can  assign  the  observed  O  to  either  H20 
or  free  O*.  Depth  profiling  of  H  will  also  permit  us  to  judge  whether  water  is 
trapped  deep  in  the  film  during  deposition  or  whether  it  diffuses  downward  after 
post-deposition  adsorption. 

*  Program  improvement  -  we  will  continue  to  improve  our  computer  program 
for  RBS  analysis. 

Simulation  of  Thermal  and  Growth  Processes 

With  the  computer  model  vastly  improved,  a  number  of  problems  will  be 
addressed  in  the  coming  year: 

*  Two  dement  films  -  the  program  has  been  modified  to  allow  particles  of 
different  sizes  to  be  deposited. 

*  Defect  propagation  -  periodic  or  irregular  defects  can  be  introduced  at  the 
substrate  or  at  higher  levels  in  the  growing  film;  the  consequences  on  further  growth 

can  be  assessed. 

*  Potentials  -  most  important,  the  simulation  has  been  made  much  more  realistic 
by  the  inclusion  of  a  potential  function  that  will  influence  the  behavior  of  adatoms 
as  they  approach  the  surface. 

Characterization  Facility 

As  in  the  past  year,  the  Measurement  Laboratory  will  serve  the  needs  of 
research  groups  in  the  Thin  Film  Center.  The  new  spectrophotometer  will  increase 
the  efficiency  of  our  measurements  and  the  accessibility  of  our  data.  Our  plans  for 
the  coming  year  include: 

*  Absolute  reflectance  accessory  -  we  intend  to  purchase  an  accessory  for  the 
Varian  2415  that  will  allow  us  to  measure  absolute  reflectance  with  high  accuracy. 

*  Electron  microscopy  -  we  will  submit  many  of  our  samples  to  analysis  in  our 
new  electron  microscope  facilities.  We  also  will  continue  our  collaboration  with  the 
Thin  Film  Group  at  Marseilles  and  the  National  Electron  Microscopy  Laboratory  at 
Arizona  State  University  in  Tempe,  Arizona. 
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Publications  citing  the  URIP  program 

J.  P.  Lehan.  J.  D.  Targove,  B.  G.  Bovard,  L.  J.  Lingg,  M.  J.  Messerly.  and  H. 

A.  Macleod,  "Intermittent  Ion  Bombardment  of  Titanium  Dioxide  Thin  Films,"  Appl. 
Opt.,  to  be  published.  (Attached  as  Appendix  A.) 

J.  D.  Targove,  "The  Ion-Assisted  Deposition  of  Optical  Thin  Films,"  Ph.D. 
Dissertation,  University  of  Arizona,  1987.  (Available  on  request.) 

Related  Publications 

Papers: 

B.  J.  Bartholomeusz.  K-H.  Muller,  and  M.  R.  Jacobson.  "Computer  Simulation 
of  the  Nucleation  and  Growth  of  Optical  Coatings,"  Proc.  SPIE  821,  to  be  published. 

B.  G.  Bovard  and  H.  A.  Macleod.  "Non-Linear  Behavior  of  Optical  Coatings 
Subjected  to  Intense  Laser  Irradiation,"  submitted  to  J.  Mod.  Opt. 

R.  Sargent  and  H.  A.  Macleod,  "Computer  Simulation  of  the  Growth  of  Two- 
Element  Films,"  Proc.  SPIE  821.  to  be  published. 

R.  Sargent  and  H.  A.  Macleod,  "Computer  Simulation  of  Substrate  Defect 
Propagation  in  Thin  Films,"  Proc.  SPIE  821,  to  be  published. 

J.  D.  Targove,  J.  P.  Lehan,  L.  J.  Lingg,  H.  A.  Macleod.  J.  A.  Leavitt,  and  L. 

C.  McIntyre.  Jr.,  Appl.  Opt.  26.  3733  (1987). 

Presentations: 

B.  G.  Bovard  and  H.  A.  Macleod,  "Thermo-Optical  Modeling  of  Multilayer 
Optical  Thin  Films  Subjected  to  Laser  Irradiation,"  at  the  1987  Annual  Meeting  of 
the  Optical  Society  of  America,  Rochester,  NY,  October  19-23. 

C-K.  Hwangbo,  L.  J.  Lingg,  J.  P.  Lehan,  J.  D.  Targove,  M.  R.  Jacobson,  H.  A. 
Macleod,  and  J.  L.  Makous,  "Ion-Assisted  Deposition  of  Thermally  Evaporated  Ag 
Films,"  at  the  1987  Annual  Meeting  of  the  Optical  Society  of  America.  Rochester, 
NY,  October  19-23. 

R.  B.  Sargent  and  H.  A.  Macleod,  "Thin  Film  Growth  Modeling,"  at  the  1987 
Annual  Meeting  of  the  Optical  Society  of  America.  Rochester.  NY,  October  19-23. 


J.  D.  Tar  gov  e,  M.  J.  Messerly,  and  H.  A.  Macleod,  "Ion-Assisted  Deposition  of 
A1,0,  Films,"  at  the  1987  Annual  Meeting  of  the  Optical  Society  of  America, 
Rochester,  NY.  October  19-23. 

F.  Varnier,  C.  Boulestix.  J.  D.  Targove,  L.  J.  Lingg.  B.  G.  Bovard,  and  H.  A. 
Macleod,  "Influence  of  Ion-Assisted  Deposition  on  Structure  and  Surface  Roughness 
of  Aluminum  Oxide,"  at  the  1987  Annual  Meeting  of  the  Optical  Society  of  America, 
Rochester.  NY,  October  19-23. 

Dissertation: 

M.  J.  Messerly.  "Ion-Beam  Analysis  of  Optical  Coatings,"  Ph.D.  Dissertation. 
University  of  Arizona,  1987. 

Personnel 


Faculty: 


H.  A.  Macleod 
M.  R.  Jacobson 
B.  G.  Bovard 


R.  H.  Potoff 

S.  Rana 
A.  McKew 

Graduate  Students: 
L.  J.  Lingg 
C-K.  Hwangbo 
J.  P.  Lehan 


Collaborations 

J.  A.  Leavitt  and  L.  C.  McIntyre,  Jr.,  Physics  Department,  University  of 
Arizona;  RBS  analysis  of  samples. 

F.  Varnier  and  C.  Boulestix,  University  of  Aix-Marseilles,  Marseilles,  France; 
electron  microscopic  investigations  of  surface  topography. 
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TASK  2B 


STRUCTURE  MODIFICATION  BY  ION-ASSISTED  DEPOSITION 
(ION-BEAM  ANALYSIS  OF  THIN  FILMS) 

Principal  Investigator:  J.  A.  Leavitt 

Project  Goals 

i 

A  principal  goal  of  this  program  is  to  provide  accurate  characterization  using 
MeV  ion-beam  analysis  (IBA)  of  thin  films  made  by  other  Thin  Film  Center 
scientists.  We  use  *He+  analysis  ion  beams,  in  the  1.5  MeV  to  5.0  MeV  energy 
!  range,  from  our  6  MV  Van  de  Graaff  accelerator  to  determine  film  thickness, 

stoichiometry,  concentration  profiles  and  crystalline  structure  by  Rutherford 
backscattering  (RBS),  high-energy  backscattering.  and  channeling. 

,  Additional  goals:  development  of  apparatus  for  film  characterization  by  nuclear 

reaction  analysis  (NRA)  of  hydrogen  in  thin  films  (using  1SN  analysis  beams), 
development  of  computer  programs  for  data  analysis,  and  investigation  use  of  heavier- 
than-4He  analysis  beams  for  RBS. 

i 

Work  Accomplished 

Film  Characterization 

» 

1.  During  the  period  9/9/86  to  10/29/87,  we  provided  complete  RBS  analyses  of  at 
least  173  thin  films  provided  by  scientists  associated  with  the  Thin  Film  Center 
(TFQ.  In  particular,  we  analyzed: 

I 

a.  139  films  for  Dr.  Macleod's  group.  Most  of  these  films  were  fluorides  of 
Nd,  Al,  Ce,  Eu,  Gd,  La,  Tb,  Ho,  Yb,  Er,  and  Sm  measured  in  connection 
with  the  use  of  ion-assisted  deposition  for  modification  of  film  structure; 
►  other  films  analyzed  were  AIN,  A1NO.  NaAlFO,  and  Ag. 


b.  12  samples  of  silica  polished  with  CeO,  grinding  compounds  for 
Oibsoo/Parks;  this  investigation  apparently  produced  no  significant 
information  concerning  surface  roughness. 

c.  17  films  with  Au  in  A1,0,  for  Suits/Gibson. 

d.  3  films  of  SiG  1  of  A1  with  N  implanted,  and  1  of  Ag  for  Salik. 

e.  2  Nb/Cu  superlattices  for  Bell/Stegeman;  we  provided  stoichiometries  of  films 
presumably  used  for  the  Brillouin  scattering  studies. 

A  portion  of  the  collaborative  work  with  the  Macleod  group  is  described  in  the 
first-listed  publication  (manuscript  attached  as  Appendix  B). 

2.  During  the  period  2/9/87  to  10/29/87,  we  provided  87  additional  complete 
hackscattering  analyses  of  films  provided  by  C.  Falco's  group  in  Physics  and  D. 
Uhlmann's  group  in  Materials  Science,  in  connection  with  their  attempts  to 
develop  superconducting  films  of  the  new  high-temperature  superconducting 
materials.  Our  analyses  provided  accurate  film  stoichiometries  and  thicknesses, 
as  well  as  information  on  film-substrate  interactions.  Falco  and  J.  Makous 
succeeded  in  making  thin  superconducting  films  of  Y,Ba1CU,07  with  a  transition 
temperature  of  75  1C  (half-height)  during  August  1987.  thus  becoming  the  third 
university  group  in  the  country  (after  Stanford  and  Cornell)  to  produce  high  Tc 
thin  films.  Our  contribution  to  this  work  was  described  in  a  poster  session.  Oct. 
1-2,  1987  at  the  Santa  Fe  meeting  (second-listed  publication,  copy  attached  as 
Appendix  Q.  Our  contribution  is  also  acknowledged  in  a  joint  publication  of 
the  Falco-Uhlmann  groups  (third-listed  publication). 

Nuclear  Reaction  Analysis  for  Hydrogen  Depth  Profiling 

Use  of  NRA  for  *H  depth  profiling  with  ISN  analysis  beams  requires  production 
of  substantial  beams  of  doubly  ionized  **N++  in  order  to  cover  the  desired 
6  MeV  -  10  MeV  energy  region  with  our  6  MV  accelerator.  We  tried  both  ,5N2  and 
,5NH,  in  our  rf  ion  source  and  found  that  ,5N++  ions  are  not  produced  in  sufficient 
quantity  for  practical  hydrogen  depth  profiling.  We  concluded  that  our  rf  ion  source 
should  be  replaced  by  a  Penning  source  (which  produces  a  larger  fraction  of  highly 
charged  ions).  We  acquired  $10,715  from  IBM-Tucson  and  bought  a  Penning  source 


and  extraction  electrode  from  General  Ionex;  at  the  time  we  did  not  have  the  funds 
for  the  recommended  einzel  lens  and  velocity  selector  to  accompany  the  source. 

We  received  the  Penning  source  and  electrode  on  S/21/87.  We  built  a  60  Hz 
bench-test  facility  and  used  it  to  obtain  data  to  determine  the  best  source  operating 
conditions,  as  well  as  whether  the  additional  components  are  necessary  to  convert  our 
accelerator  into  a  practical  NRA  device.  As  a  result  of  the  tests,  we  have  concluded 
that  we  do  need  the  additional  components.  We  converted  $10,295  in  capital  funds 
in  our  section  of  the  TFC  contract  for  purchase  of  the  velocity  selector  and  have 
acquired  $6,910  from  the  Optical  Data  Storage  Center  for  purchase  of  the  einzel  lens. 
These  items  were  ordered  10/1/87.  We  anticipate  installation  in  December  1987. 

Other 

1.  We  have  continued  work  with  heavier-than-‘He  analysis  beams  to  improve  mass 
and  depth  resolution  in  RBS.  Preliminary  analysis  of  recent  data  obtained  with 
C.  At,  and  Kr  analysis  beams  indicates  no  major  improvement  in  mass  and  depth 
resolution  will  result  from  using  the  heavier  beams. 

2.  We  have  continued  to  develop  least-squares  fitting  programs  for  accurate 
deconvolution  of  overlapping  backscattering  peaks.  We  have  checked  the 
reliability  of  the  results  produced  by  the  deconvolution  by  comparing  with  results 
obtained  from  data  run  at  energies  sufficient  to  produce  separated  peaks.  See 
the  fourth-listed  publication  below  (manuscript  attached  as  Appendix  D). 

3.  On  10/23/87,  a  research  briefing  by  J.  Leavitt  describing  our  work,  entitled 
MMeV  Ion  Beam  Analysis  of  Thin  Films,"  was  broadcast  live,  nationwide,  by 
satellite.  This  briefing  contained  short  discussions  of  analysis  techniques,  with 
illustrative  examples  from  the  fields  of  optical  coatings,  high  Tc  superconducting 
films,  optical  data  storage.  A  video  tape  of  this  one-hour  briefing  may  be 
obtained  from  Microcampus,  University  of  Arizona. 

Plans  for  Coming  Year 

A.  We  plan  to  continue  film  characterization  for  our  TFC  collaborators  (and 
others).  We  expect  the  sample  load  to  increase  considerably  during  the  this  next 
year  as  film-making  plants  commence  operation. 
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B.  We  will  exert  a  major  effort  starting  in  December  1987  to  bring  the  ‘H 
depth  profiling  capability  using  ISN  beams  on  line.  We  expect  to  require 
approximately  two  months  to  get  the  Penning  source  installed  and  operating  smoothly. 

C.  We  plan  to  measure  non-Rutherford  cross  sections  for  the  scattering  of  He 
from  light  elements  in  the  2  MeV  to  10  MeV  region. 

D.  We  expect  to  continue  development  of  our  computer  fitting  programs  for 
improved  data  analysis. 

Publications 

J.  D.  Tar  gov  e,  L.  J.  Lingg.  B.  G.  Bovard,  J.  P.  Lehan.  H.  A.  Macleod,  J.  A. 
Leavitt,  and  L.  C.  McIntyre,  Jr.,  "Preparation  of  Aluminum  Nitride  and  Oxynitride 
Thin  Films  by  Ion-Assisted  Deposition,"  to  appear  in  MRS  Symposia  Proceedings  - 
Materials  Modification  and  Growth  Using  Ion  Beams  -  Symposium  C.  1987  Spring 
Meeting  (copy  attached.  Appendix  B). 

Poster  session  paper: 

J.  A.  Leavitt  and  L.  C.  McIntyre,  Jr.,  "Stoichiometry  of  High  Tc  Thin  Filins  by 
High  Energy  Backscattering,"  at  New  Mexico  Section  MRS  Conference  on 
Developments  in  High-Temperature  Supeconducting  Materials,  Oct.  1-2.  1987,  Santa 
Fe,  NM  (copy  attached.  Appendix  Q. 
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Uhlmann.  and  D.  R.  Uhlmann,  "Superconducting  and  Structural  Properties  of 
Sputtered  Thin  Films  of  Y  Ba,Cus07_x,"  submitted  to  Appl.  Phys.  Lett. 

L.  C.  McIntyre,  Jr.,  M.  D.  Ashbaugh,  and  J.  A.  Leavitt,  "Limits  on  Accuracy  of 
Stoichiometry  Determined  by  Rutherford  Backscattering  Using  Computer  Peak  Fitting," 
to  appear  in  MRS  Symposia  Proceedings  -  Materials  Modification  and  Growth  Using 
Ion  -  Symposium  C,  1987  Spring  Meeting  (copy  attached.  Appendix  D). 
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Video  Tape: 


J.  A.  Leavitt,  "MeV  Ion  Beam  Analysis  of  Thin  Films,"  National  Technological 
University/University  of  Arizona  Microcampus  Research  Seminar,  live  satellite 
broadcast /later  video  tape  distribution.  60  minutes,  Oct.  23,  1987  (video  tape  may  be 
obtained  from  UA  Microcampus). 

Personnel 


J.  A.  Leavitt 
P.  Stoss 

M.  D.  Ashbaugh 
Collaborations 


All  our  work  is  collaborative.  See  collaborators  listed  above. 


TASK  3 


A  NOVEL  TECHNIQUE  FOR  QUANTIFYING 
THE  MECHANICAL  PROPERTIES  OF  THIN  FILMS 

Principal  Investigators:  G.  I.  Stegeman  and  C.  T.  Seaton 

Project  Goals 

The  principal  objective  of  this  task  is  to  develop  Brillouin  scattering  as  a 
technique  for  measuring  the  elastic  properties  of  thin  films  and  interfaces.  The  work 
performed  can  be  subdivided  into  two  areas: 

a)  volume  elastic  properties  of  films  and  their  modification; 

b)  interface-induced  modifications  of  film  properties. 

Summary 

We  developed  the  experimental  tools  and  analysis  techniques  for  determining 
elastic  constants  of  films  by  Brillouin  spectroscopy.  Stonley  and  Love  waves  were 
observed  for  the  first  time  with  Brillouin  scattering.  The  Stonley  waves  were  used 
to  compare  the  mechanical  properties  of  a  tantalum  film  at  its  air  and  substrate 
interfaces.  The  Love  waves  were  used  to  measure  the  c,3  elastic  constant  for  films 
for  the  first  time.  In  addition,  new  longitudinally  polarized  thin-film  guided  waves 
were  observed  and  used  to  determine  (in  some  cases)  the  cn  elastic  constants  with 
much  higher  precison  than  previously  possible.  Finally,  modifications  in  the  elastic 
constants  upon  polymerization  of  Langmuir-Blodgett  films  are  being  measured. 

We  have  concentrated  on  modeling  and  measuring  interface  effects  using 
Brillouin  scattering.  As  a  precursor  to  the  experiments,  we  modeled  theoretically  the 
adhesion  of  thin  films,  and  found  acoustic  dispersion  characteristics  sensitive  to  film 
adhesion.  We  measured  films  provided  by  IBM-Tucson  to  test  the  sensitivity  of 
Brillouin  scattering  to  interface  stresses,  without  success.  Similar  investigations  of 
metallic  superlattices  with  large  "interface  densities"  were  successful.  As  a 
continuation  of  the  interests  of  B.  Hillebrands.  we  carried  out  Brillouin  scattering 


experiments  on  spin  waves  in  magnetic  films  and  superlattices.  Interface  magnetism, 
which  reveals  important  information  about  the  interdiffusion,  electronic  hybridization 
and/or  polarization  at  interfaces,  as  well  as  about  interface  anisotropies  was 
investigated. 

Work  Accomplished 

In  the  following  sections  we  report  the  recent  results  of  this  project  and  outline 
future  investigations. 

Determination  of  Film  Elastic  Constants 

We  developed  computer  programs  fen-  deducing  elastic  constants  of  thin  films 
from  experimentally  obtained  surface-phonon  dispersion  curves.  The  fitting 
procedure  requires  that  the  transverse  sound  velocity  of  the  film  be  significantly 
smaller  than  that  of  the  substrate,  which  normally  can  be  achieved  by  a  suitable 
choice  of  substrate  material.  The  existing  programs  also  require  that  the  symmetries 
of  the  film  and  the  substrate  are  at  least  tetragonal.  Furthermore,  at  minimum,  either 
the  film  or  the  substrate  must  be  opaque.  The  best  precision  can  be  achieved  for  a 
sequence  of  films  with  thicknesses  ranging  from  100  A  to  5000  A.  However,  in 
some  special  cases,  the  elastic  constants  could  be  obtained  for  single  films.  Elastic 
constants  have  been  determined  for  ZnSe(00I)  films  on  GaAstOOl).1-*  CdA-  and  ODF- 
Langmuir-Blodgett  films.1*4  Ta  films.1  Cu/Nb.**T  and  Fe/Pd  films.*-*  In  many  cases 
the  total  set  of  elastic  constants  was  determined. 

Evidence  for  the  Existence  of  Guided  Longitudinal 
Acoustic  Phonons  in  ZnSe  on  GaAs 

We  found  theoretical  as  well  as  experimental  evidence  for  a  new  class  of  long- 
wavelength  acoustic  film  excitations  which  we  identified  as  longitudinal  guided 
acoustic  phonons.1  Despite  their  coupling  to  propagating  transverse  film  modes  at  the 
film  boundaries,  these  new  modes  have  lifetimes  comparable  to  the  Rayleigh  mode. 
They  are  the  longitudinal  counterpart  of  the  well-known  Lamb  or  Sezawa  modes, 
with  analogous  existence  criteria.  The  characteristic  necessary  for  observation  is  a 
transparent  film  at  least  1  fim  thick,  whose  longitudinal  sound  velocities  are  smaller 
than  those  of  the  substrate.  It  is  most  important  that  these  films  permit  a  highly 
accurate  determination  of  c,,  for  the  film,  a  parameter  which  is  otherwise  difficult  to 


evaluate.  Since  the  displacement  field  of  the  modes  is  nearly  completely  polarized 
parallel  to  the  film,  the  mode  can  also  be  used  to  monitor  effects  which  modulate 
only  dilatational  waves. 

Elastic  Properties  of  Metallic  Films  and  Interfaces 

We  studied  Brillouin  scattering  from  a  number  of  metallic  film  systems, 
including  tantalum  films  and  Cu/Nb  and  Fe/Pd  superlattices.  Superlattices  are  very 
well  suited  for  studying  interface-related  effects  because  of  their  large  number  of 
interfaces.  We  measured  the  elastic  constants  for  Nb/Cu  and  Fe/Pd  superlattices.7*8 
The  elastic  constants  of  Cu/Nb  and  Fe/Pd  superlattices  exhibit  a  large  softening  of 
cu  for  small  modulation  wavelengths.  The  results  are  indicative  of  the  importance  of 
the  band-structure  properties  on  the  elastic  properties.  A  hybridization  of  the 
electronic  states  at  the  interface  and  a  corresponding  modification  of  the  density  of 
states  at  the  Fermi  level  in  the  interface  region  (hybridization  length)  presumably 
accounts  for  the  observed  behavior.  A  more  quantitative  understanding  of  the 
observed  phenomena  is  strongly  demanded. 

The  cu  elastic  constant  of  the  Cu/Nb  superlattice  was  obtained  by  measuring 
the  Love  acoustic  modes  supported  by  the  film  for  the  first  time  with  Brillouin 
scattering.*  The  salient  question  is  whether  there  is  also  a  softening  of  the  c12  elastic 
constant  with  modulation  wavelength.  None  was  measured  experimentally  to  within 
an  accuracy  of  10%. 

Molybdenum  films  on  sapphire  substrates  satisfy  the  conditions  for  the  existence 
of  Stonley  waves.1  These  modes  are  highly  localized  at  the  film-substrate  interface, 
in  contrast  to  Rayleigh  waves  which  are  localized  at  the  film-air  interface. 
Measurements  of  the  dispersion  curves  for  these  two  modes  on  carefully  prepared 
films  showed  that  the  tantalum  film  elastic  constants  at  both  film  boundaries  are 
essentially  identical. 

Elastic  Properties  of  CdA  and  ODF  Langmuir- Blodgett  Films  and  their 
Modifications  upon  Film  Polymerization 

We  studied  the  elastic  properties  of  CdA-,  saturated  ODF-  and  unpolymerized 
and  polymerized  w-ODF-  and  ODF:ODM-Langmuir-Blodgett  films.8*4  Because  of  the 
laminar  film  structure  consisting  of  parallel  aligned,  long  molecules,  which  are 


perpendicular  to  the  interface  and  are  tightly  bound  to  each  other,  the  films  are 
elastically  highly  anisotropic  and  show  generally  an  unusually  small  value  of  the 
transverse  sound  velocity,  i.e.,  cu.  For  CdA-Langmuir-Blodgett  films,  all  of  the 
elastic  constants  could  be  determined.4  The  first  preliminary  results  for  the  ODF 
systems  indicate  a  very  puzzling  behavior  of  the  elastic  constants  upon 
polymerization:  whereas  the  ODF:ODM  films  show  a  small  stiffening  of  cM  upon 
polymerization,  the  u>-ODF  samples  show  a  decrease  in  the  Rayleigh  sound  velocity 
and  a  slight  stiffening  of  the  Sezawa  modes.  This  also  indicates  large  modifications 
of  the  longitudinal  sound  velocity,  which  are  not  easily  measured  in  a  direct  way. 
Future  investigations  should  permit  a  more  detailed  understanding  of  the  elastic 
properties  of  these  films. 

Modeling  Adhesion  by  Infinitely  Thin  Buffer  Layers: 

New  Long-Wavelength  Acoustic  Film  Excitations 

We  used  computer  simulations  to  study  the  influence  of  adhesion  on  the 
dispersion  curves  of  long- wavelength  film  excitations.  We  modeled  the  interface  as  a 
very  thin  buffer  layer,  characterized  by  a  vanishing  transverse  sound  velocity.  In 
regions  where  the  dispersion  curves  show  a  nearly  zero  slope  without  a  buffer  layer, 
a  new  mode  splits  off  in  the  presence  of  the  buffer  layer,  which  is  reminiscent  of 

the  full  symmetric  mode  of  a  free-standing  film.  Future  investigations  will  show 

whether  the  predicted  new  modes  can  be  found  experimentally,  and  whether  they 
will  allow  us  to  characterize  adhesion. 

Characterization  of  Interface  Stresses 

We  attempted  to  measure  the  influence  of  interface  stresses  on  elastic  constants. 
For  ZnSe  films  on  GaAs  it  is  known  that  interface  stresses  exist  in  a  region  up  to 
1  fim  from  the  interface,  before  they  are  relaxed  by  the  formation  of  dislocations. 
Since  stress  fields  do  not  enter  the  acoustic  equations  of  motion  directly,  the  only 
means  of  detecting  them  with  Brillouin  spectroscopy  is  by  a  renormalization  of  the 

elastic  constants  due  to  higher-order  terms.  We  could  not  find  any  thickness 

dependence  of  the  elastic  constants  within  the  experimental  accuracy  of  about  5%  for 
this  system.1 


Interfact  Magnetism 


In  studies  of  the  structural  and  electronic  properties  of  interfaces,  investigation 
of  interface  magnetism  can  shed  light  on  many  questions.  In  particular,  the 
sensitivity  of  magnetic  moments  on  impurities  (interdiffusion)  and  on  modifications  of 
the  electronic  band  structure  (hybridization,  charge  transfer),  as  well  as  film  growth- 
induced  magnetic  interface  and  volume  anisotropies  are  excellent  tools  for  the 
characterization  of  film  properties.9*10  Studying  magnetism  by  measuring  spin-wave 
spectra  in  a  Brillouin  scattering  experiment  allows  a  quantitative  determination  of  the 
saturation  magnetization,  the  g-f actor,  the  exchange  constant,  volume  and  interface 
anisotropies,  the  interface  interdiffusion  range,  and  possible  magnetic  spacer 
polarizations.  We  observed  and  analyzed  collective  magnetostatic  spin-wave 
excitations  in  magnetic  Fe/Pd  and  Fe/W  superlattices.  For  Fe/Pd.  we  found  strong 
evidence  for  a  magnetic  polarization  of  the  Pd  spacer  layers,  as  well  as  for  a  small 
negative  out-of-plane  interface  anisotropy  constant.  The  reference  system  Fe/W 
shows  no  pronounced  anomalies. 

This  portion  of  the  project  was  not  part  of  the  original  program  and  represents 
an  extension  of  previous  interests  for  Dr.  Burkard  Hillebrands. 

Plans  for  the  Coming  Year 

The  following  investigations  are  tentatively  scheduled  for  the  coming  year: 

1.  Characterization  of  the  elastic  properties  of  ZnS  films. 

2.  Elastic  properties  of  Langmuir-Blodgett  films  upon  polymerization. 

3.  The  adhesion  problem. 

4.  Characterization  of  surface  roughness  by  measuring  surface  phonons. 

5.  Characterization  of  magnetic  interface  properties. 

Topic  1  is  an  application  of  the  described  new  method  for  determination  of 
transparent  film  elastic  constants.  Films  grown  in  Task  1  by  various  methods  will  be 
characterized  by  Brillouin  spectroscopy.  Topic  2  is  a  continuation  of  the  current 
investigations.  In  Topic  3  we  will  test  our  theoretical  predictions.  It  is  highly 
possible  that  the  Langmuir-Blodgett  films  will  provide  good  testing  systems.  Topic  4 
is  intended  to  check  theoretical  predictions  about  roughness-induced  frequency 
renormalizations  of  the  Rayleigh  mode.  The  method  will  allow  characterization  of 


roughness  of  highly  corrugated  surfaces  on  a  length  scale  of  surface-phonon 
wavelengths.  In  Topic  S  we  will  investigate  in  more  detail  the  Fe/Pd  interface  and 
attempt  to  correlate  the  observed  elastic  and  magnetic  anomalies. 

Further  additional  possibilities  exist.  In  particular: 

1.  A  search  for  a  second  surface  mode  on  viscoelastic  media. 

2.  Correlation  of  elastic  film  properties  with  material  hardening  and  working,  in 
collaboration  with  IBM-San  Jose. 

Personnel 

G.  I.  Stegeman 

B.  Hillebrands 
S.  Lee 

Collaborations 

U.  J.  Gibson.  UA  -  adhesion  properties  of  thin  films. 

C.  M.  Falco,  UA  -  preparation  of  superlattice  samples. 

J.  E.  Potts.  3M  Company.  St.  Paul.  MN  -  preparation  of  ZnSe  samples. 

J.  Swalen,  IBM.  San  Jose,  CA  -  preparation  of  Langmuir-Blodgett  films. 
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TASK  4 


CONTINUUM  MODEL  OP  FILM  GROWTH 
Principal  Investigator:  Seth  Lichter 
Project  Goal 

Develop  a  theoretical  continuum  approach  to  modeling  the  microstructure  of  thin 
films. 

Summary 

Most  thin-film  studies  use  a  Monte-Carlo  simulation  of  the  deposition  process. 
These  models  omit  much  of  the  relevant  thermodynamics  and  dynamics  which  govern 
adatom  mobility.  This  omission  is  necessitated  by  the  complexity  and  size  of  the 
simulation.  Those  few  simulations  which  account  for  particle-particle  interactions  are 
limited  to  extremely  small  grids. 

Additionally,  these  simulations  do  not  cope  with  another  fundamental  aspect  of 
thin-film  physics:  the  deposited  film  must  satisfy  some  constitutive  constraints  as 
well  as  constraints  on  stress  and  strain. 

It  is  our  intention  to  pursue  an  alternative  to  the  usual  numerical  computation. 
Rather  than  starting  with  numerical  simulation,  we  will  begin  with  a  simulation  of 
the  governing  physical  processes.  We  will  preserve  the  physics  as  much  as  possible, 
in  the  belief  that  a  properly  formulated  model  will  yield  correct  results. 

The  one  assumption  common  to  all  our  work  is  that  the  continuum  hypothesis 
can  be  applied.  We  recently  applied  this  approach  to  a  model  of  adatom  mobility 
which  gave  an  analytical  basis  for  the  tangent  rule.1  While  that  initial  work  focused 
on  the  deposition  process,  this  complementary  effort  will  focus  on  the  equations 
describing  the  deformation  of  the  film  itself. 


Work  Accomplished 


In  our  work,  we  have  ignored  the  deposition  process  entirely,  and  have  solved 
for  possible  equilibrium  structures  in  the  films.  Ultimately  we  hope  to  use  this 
model  to  predict  the  abrupt  structural  transitions  catalogued  by  Movchan  and 
Demchishin.*  Our  attempts  in  this  area  begin  with  analyses  of  one-dimensional 
deformable  lattices.  This  work  is  in  its  initial  stages. 

Other  work  has  centered  on  the  motion  of  grain  boundaries  and,  in  particular, 
those  conditions  under  which  the  grain  boundaries  could  organize  themselves  into 
columnar  microstructures.  This  work  is  still  in  the  tentative  stage,  but  we  show 
below  how  a  simple  model,  based  on  the  usual  equations  of  mechanical  equilibrium 
in  an  elastic  thin  film,  can  be  used  to  predict  column  size  as  a  function  of  tensile 
stress  in  the  film. 

The  three-dimensional  displacements  for  an  elastic  body  are  governed  by  the 
equation 


d*Ui 

/>  a?-  -  Fi  +  /*uMi  +  (^  +  m)  Ujjj  .  (i) 

subject  to  the  boundary  conditions 

Uj  -  fj  (X;,  t)  .  (2) 

where  f;  is  continuous,  u;  is  the  displacement,  Fj  is  the  body  force,  p  is  the  shear 
modulus,  p  is  the  density,  and  the  subscripts  following  a  comma  denote  derivatives. 
We  assume  the  flux  of  grain  boundaries  B-,  is  proportional  to  the  body  force 

B;  -  D  F;  .  (3) 

where  D  is  a  proportionality  constant.  We  further  assume  that  the  divergence  of  the 
flux  of  grain  boundaries  is  compensated  for  by  a  creation  C  and  an  annihilation  A 
of  boundaries  such  that 


-D  V*  U  -  C  -  A  . 
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(4) 


Subtracting  the  elastic  solution,  restricting  the  force  to  lie  in  a  thin  plane  in  the 
xs  direction,  and  after  some  calculation  we  find 


[p  J?~ 

M  v*  j  v>  u,*  -  -v,1  uy  . 

(5a) 

-D  V*  U'  -  fi  uM*  . 

(5b) 

Uj'  -  fj'  on  the  boundaries. 

(5c) 

where  the  primes  denote  the  total  displacement  minus  the  elastic  displacement.  A 
particular  solution  can  be  found  by  suppressing  the  time  dependence.  One 
approximate  solution  is  then 


w  -  w0 


where  fi  and  w0  are  constants  and  a  is  the  characteristic  horizontal  wave  number. 
As  w  has  the  same  dependence  on  the  applied  stress  as  does  the  displacement  due  to 
motion  of  the  grain  boundaries,  we  find 


T“  *  a(^]‘/‘  *  (7) 

where  ru  is  the  stress.  So  the  solution  reveals  a  periodic  array  of  grain  boundaries. 
Such  an  array  may  delineate  a  row  of  the  columnar  microstructure.  Furthermore,  in 
this  limit,  the  size  of  the  columns  is  shown  to  be  inversely  proportional  to  the  stress 
in  the  film. 
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Plans  for  the  Coming  Year 


This  year  we  hope  to  show  that  the  possible  equilibrium  solutions  are  dependent  on 
the  conditions  in  the  film.  In  the  above  example  the  controlling  parameter  is  film 
stress.  The  stress  can  be  related  to  substrate  temperature  using  energy 
considerations.  In  this  manner  we  hope  to  establish  the  dependence  of  columnar 
microstructure  on  substrate  temperature. 
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S.  Lichter 
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TASK  6 


SURFACES  AND  POLISHING  STUDIES 
Principal  Investigators:  R.  E.  Sumner  and  R.  E.  Parks 
Project  Goals 

An  important  aspect  of  the  development  of  thin  films  is  the  manufacture  and 
characterization  of  known  substrates.  The  objectives  of  the  research  program  in 
substrate  production,  conducted  at  the  Optical  Shop,  are: 

1.  Define  and  quantify  the  different  surface  defects  that  may  be  present  in  a 
substrate  surface. 

2.  Determine  the  surface  quality  achieved  by  different  polishing  procedures  in  a 
number  of  materials  (glasses,  crystals,  and  metals). 

3.  Define  efficient  polishing  procedures  to  produce  surfaces  with  a  given  quality, 
with  special  emphasis  on  the  production  of  low-scatter,  smooth  surfaces. 

Work  Accomplished 

Goal  1. 

Initial  library  research  has  been  conducted. 

An  inspection  optical  microscope  to  identify  visible  surface  defects  has  been 
installed. 

A  non-contact  optical  profilometer  to  measure  surface  roughness  has  been 
installed. 

A  scatterometer  is  under  construction  (see  Fig.  I). 

Goal  2. 

As  a  preliminary  step  in  designing  the  grinding  and  polishing  procedures,  some 
simple  experiments,  that  is,  observation  of  the  grinding  and  polishing  of  quartz 
substrates,  were  conducted.  The  results  of  these  experiments  suggest  criteria  to 
consider  during  the  grinding  process,  show  how  critical  the  cleaning  of  the  samples 
may  be  when  characterizing  them,  and  also  suggest  that  a  variety  of  testing 


instruments  should  be  used  to  better  characterize  and  understand  the  polishing 
process. 


At  this  point  in  time  some  of  the  necessary  materials  have  been  ordered  and 
received,  quartz  samples  and  their  supports  have  been  prepared,  the  milling  of  the 
polishing  materials  has  been  started,  flat  cast-iron  grinding  tools  have  been  prepared, 
and  a  measurement  table  has  been  constructed  to  quantify  material  removal  during 
the  grinding  process. 

Two  blocks  of  19  quartz  samples  each  have  been  ground  and  polished  using 
two  different  procedures.  Their  characterization  is  under  way.  Two  other  blocks 
are  being  prepared. 

Notes  on  the  grinding  and  polishing  procedures  are  being  taken. 


Figure  1.  A  diagram  of  the  scatterometer. 


Plans  for  the  Coming  Year 


Several  blocks  of  quartz  samples  will  be  prepared  using  different  polishing 
procedures. 

Several  blocks  of  BK-7  glass  samples  and  crystal  samples  will  be  produced. 
Polishing  procedures  will  be  directed  by  past  experience  with  the  quartz  samples. 

The  surface  quality  will  be  characterized  for  the  blocks  produced. 

Personnel 

R.  E.  Summer 
R.  E.  Parks 
J.  M.  Sasian 

Collaborations 

A  number  of  glass  substrates  have  been  prepared  for  other  members  of  the 
group. 
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TASK  7 


CHARACTERIZATION  FACILITY 
Project  Goal 

Set  up  a  central  facility  for  characterization  of  thin  films,  surfaces,  and 
substrates. 

Work  Accomplished 

Major  equipment  purchases  and  hiring  of  a  full-time  technician  to  oversee  the 
facility  were  the  major  accomplishments  in  this  year.  The  transmission  electron 
microscope  (TEM)  and  the  scanning  electron  microscope  (SEM)  were  installed  in 
October,  and  the  instruments  are  operational.  Not  all  of  the  sample  preparation 
equipment  has  been  received.  The  WYKO  profilometer  has  been  operational  for 
some  time,  and  is  being  used  for  both  thickness  determinations  and  surface 
characterization.  A  series  of  seminars  on  the  equipment  capabilities  and  limitations 
is  now  underway. 

A  small  effort  was  initiated  under  this  task  to  tap  local  expertise  on  scanning 
tunneling  microscopy.  We  expect  that  a  vacuum-compatible  working  model  will  be 
completed  in  the  next  year. 

Acquisitions  for  Characterization  Facility 

1.  JEOL  2000 FX  Transmission  Electron  Microscope: 
o 

1.4  A  line  resolution, 

bright  field,  dark  field,  diffraction,  and  nano-diffraction  modes. . 

Sample  preparation  equipment  for  TEM: 
ion  thinner  for  cross-sectional  samples, 
dimpler  for  cross-sectional  samples, 
chemical  jet  thinner  for  cross-sectional  samples, 
precision  core  drill  for  cross-sectional  samples. 


2.  JEOL  JSM-U3  Scanning  Electron  Microscope: 

60  A  resolution. 

Sample  preparation  equipment  for  SEM: 
sputter  coater,  target,  etc. 

3.  WYKO  2-d  heterodyne  interferometer  (on  loan). 
Other  equipment: 

4.  Varian  CARY  2415  spectrophotometer, 

donated  to  Dr.  Macleod's  Measurements  Lab. 

5.  Scatterometer. 

built  as  part  of  Task  6. 

Personnel 

U.  J.  Gibson 

V.  Lindley 
M.  Kelley 
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Intermittent  Ion  Bombardment 
of  Titanium  Dioxide  Thin  Films 

J.  P.  Lehan,  J.  0.  Targove,  B.  0.  Bovard,  L  J.  Lingg,  M.  J.  Maw  arty, 
and  H.  A.  Mad  sod 

Titanium  dioxida  is  important  to  tha  optical  coating  industry  bacausa  of  its  high 
indax;  unfortunately.  it  is  not  aasily  deposited  as  stoichiometric  TiO*.1  Oxygen 
deficiency  in  TiO,  results  in  a  brood  absorption  bond,  cantered  about  1.05  mi n.  first 
identified  in  single  crystals  by  Cronemeyer*  and  later  observed  in  evaporated  films 
by  Hermann.*  Tha  conventional  method  of  depositing  TIO*  (which  refers  to  the 
oxides  of  11  with  stoichiometries  near  tha  desired  TiO*).  is  reactive  evaporation  onto 
a  substrata  heated  to  about  200*0  or  mors4  in  an  oxygen  atmosphere  to  promote  film 
oxidation.  Tha  elevated  temperature  hardens  the  coatings;  however,  it  increases 
absorption.1  In  addition,  high  temperatures  preclude  deposition  on  a  variety  of 
substrates,  including  plastics  and  high-precision  optical  flats.  A  low-temperature 
alternative  to  reactive  evaporation  ia  ion-assisted  deposition  (LAD),  which  is  an 
effective  method  of  improving  the  adhesion.1  stress.1  pocking  density.*  and 
environmental  stability*  of  a  wide  variety  of  costing  materials.  Martin.7  McNeil*  and 
Allen*  have  reported  on  the  LAD  of  TiO*  films.  These  studies,  however,  along  with 
most  other  work  on  the  subject,  deal  with  stationary  substrates  under  a  constant  ion 
flux.  This  would  not  be  the  situation  in  a  production  environment,  where  planetary 
rotation  ensures  uniform  coverage  of  many  substrates.  To  employ  LAD  in  such  an 
environment,  the  ion  beam  must  move  synchronously  with  the  substrates  or  the  beam 
must  extend  far  enough  to  immerse  all  of  the  substrates  in  a  constant  ion  flux. 
Alternatively,  the  substrates  can  be  rotated  through  a  stationary  smaller  ion  beam, 
making  the  bombardment  intermittent.  The  last  approach  is  the  easiest  and  requires 
the  least  modification  of  an  existing  chamber. 

J.  P.  Lehan.  B.  O.  Bovard,  L.  J.  Lingg,  and  H.  A.  Macleod  are  with  the  Optical 
Sciences  Center.  University  of  Arizona.  Tucson.  Arizona  85721.  J.  D.  Targove  and 
M.  J.  Mtsserly  are  both  formerly  with  Optical  Sciences.  J.  D.  Targove  is  presently 
with  the  Air  Force  Institute  of  Technology.  Department  of  Engineering  Physics, 
Wright-Patterson  AFB  and  M.  J.  Messerly  is  with  3M. 
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We  investigated  the  consequences  of  intermittent  ion  bombardment  on  TiO,  films 
as  the  substrate  rotated  through  the  ion  beam.  Cerac™  Ti30,  tablets  were  electron- 
beam  evaporated  in  a  cryopumped  Balzers  BAK  760  box  coaler  equipped  with  a 
3  cm  aperture  Commonwealth  Scientific™  Kaufman  hot-cathode  ion  source.  We 
added  oxygen  to  the  chamber  both  through  the  ion  gun  and  the  side  of  the  chamber 
ae  a  neutral  background  species,  to  aid  oxidation  of  the  growing  films.  The 
deposition  rate  and  Him  thickness  were  monitored  by  a  quartz  crystal  oscillator.  The 
transmittance  and  reflectance  of  the  coatings  were  measured  by  a  Cary  14 
spectrophotometer.  We  computed  the  optical  constants  from  the  tramsmittance  and 
reflectance  profiles  using  an  envelope  calculation  based  on  an  inhomogeneous,  weakly 
absorbing  film  model  developed  by  Madeod.1* 

The  chamber  parameters  during  the  deposition  were 

Substrate  temperature:  25*0,  rising  to  150*C  during  deposition 
Chamber  base  pressure:  2x10'*  mbar 
Added  gases:  O, 

Pressure:  8x10**  mbar  admitted  through  the  ion  gun 

2x1 0“*  mbar  through  the  side  of  the  chamber 
Ion  beam  voltage:  300  V 
Ion  beam  duty  cycle:  *50% 

Evaporation  rate:  4  A/s 

Ion  gun-to-substrate  distance  (minimum):  38  cm. 

The  ion-beam  center  was  aimed  at  the  point  retrieving  the  maximum  evaporant  flux 
along  the  substrate's  rotation  path,  which  was  directly  over  the  electron  beam 
crucible.  Therefore,  the  maximum  evaporant  and  ion  fluxes  coincided  at  the 
substrate.  The  ions  impinged  on  the  substrates  at  approximately  normal  incidence 
while  the  evaporant  arrived  at  about  15*  from  the  normal.  Despite  the  small  size  of 
the  ion-gun  aperture,  the  resultant  ion  beam  was  very  wide,  about  60  cm  full  width 
at  half-maximum  current  density  as  measured  along  the  substrate  rotation  path.  This 
curious  effect  is  probably  attributable  to  ion  scattering  from  the  relatively  abundant 
oxygen  in  the  chamber.  At  3x10“*  mbar  and  125°C  the  mean  free  path  of  an 


oxygen  molecule  is  13.7  cm.  less  than  our  ion  gun-to-substrate  distance  of  38  cm. 
We  later  reduced  our  ion  gun-to-substrate  distance  to  30  cm  to  reduce  the  scattering. 
The  data  were  taken  with  300  eV  oxygen  ions,  as  500  eV  ions  damaged  the  films 
even  at  modest  ion  current  densities. 
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The  number  of  ions  per  arriving  molecule.  7,  has  been  postulated  as  the 
physically  significant  parameter  for  IAD.*  With  rotation  of  the  substrates,  however, 
7  is  a  function  of  the  substrate  location.  We  thus  propose  that  the  time-averaged 

ratio.  7.  be  substituted.  To  calculate  7  we  measured  the  ion-beam  profile  along  the 
substrate  rotation  path  with  a  Faraday  cup  biased  at  -45  V  to  repel  the  neutralizing 
electrons.  The  measured  ion  current  density  was  then  fitted  to  a  function;  the  best 
fit  to  the  measured  data  was  found  to  be  Lorentzian.  The  relationship  between  the 

current  density  and  the  measured  evaporation  rate  allowed  us  to  calculate  7.  Since 
our  rotation  was  along  a  circular  path  and  the  rotation  speed  was  constant,  the  ion 
flux  was  periodic  and  so  averaged  over  a  single  period.  Then 
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We  assume  that  the  radius  is  constant  for  the  circular  rotation  and  that  t<0)  is  a  best 
fit  This  parameter  correlated  well  with  the  observed  changes  in  film  properties 
whereas  the  peak  in  7(6)  did  not. 

Because  of  the  strong  refractive-index  inhomogeneity  observed  in  some  of  our 

TiO,  films,  we  report  both  an  average  index,  n.  and  the  degree  of  inhomogeneity. 

An/n.  Figure  1  shows  n  at  550  nm  as  a  function  of  7.  Notice  that,  at  an  ion  gun- 

to-substrate  distance  of  38  cm,  n  increases  until  it  reaches  a  maximum  at  7  at  0.5. 
At  an  ion  gun-to-substrate  distance  of  30  cm.  the  slope  of  the  curve  is  greater.  This 
may  be  attributable  to  the  reduced  scattering  of  ions  by  the  oxygen  backfill.  In 

Fig.  2.  we  plot  the  extinction  coefficient,  k.  as  a  function  of  7.  Here  Ic  is  minimum 

in  roughly  the  same  7  region  as  the  maximum  in  n.  Note  that  the  film  with  7  a  0.4 

in  Fig.  2  corresponds  to  an  extinction  coefficient  value  below  the  resolution  of  our 
spectrophotometric  measurements  and  so  was  estimated  for  the  plot. 

Figure  3  plots  An/n  at  550  nm  as  a  function  of  7.  The  inhomogeneity  is  taken 
as  negative  if  the  index  decreases  away  from  the  substrate.  The  critical  point 

demonstrated  here  is  that  the  minimum  degree  of  inhomogeneity  occurs  at  7  a  0.4. 

the  same  value  that  yielded  a  maximum  in  refractive  index  and  a  minimum  in 

extinction  coefficient.  We  also  found  that  the  dispersion  of  the  degree  of 
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inhomogeneity.  An/n  vs  X.  is  altered  by  intermittent  IAD.  At  7  -  0.4S  we  observed 
an  approximately  linear  dispersion  curve  for  the  degree  of  inhomogeneity  in  the 
visible  region  of  the  spectrum  described  by  the  equation 

An/n  -  -0.0001  IX  (nm)  ♦  0.0704. 

We  gain  insight  into  the  underlying  physical  processes  if  we  examine  the 
vacuum-to-air  shift  (the  shift  in  the  spectral  profile  as  the  evacuated  chamber  is 
brought  up  to  atmospheric  pressure)  at  a  wavelength  at  which  the  optical  thickness  of 
the  coating  is  an  integral  number  of  half  wavelengths.  If  we  assume  that  this  shift 
is  caused  by  water  adsorption  into  the  film  pores,7  this  provides  a  rough  measure  of 

their  porosity.  Figure  4  plots  the  vacuum-to-air  shift  against  7.  For  7  below  0.40. 
intermittent  ion-assisted  deposition  does  not  densify  the  films  significantly,  but  above 

this  level  some  densiflcation  of  the  films  occurs.  Above  7  a  0.45.  however, 
absorption  rises  steeply  and  index  decreases. 

Other  analysis  techniques  included  x-ray  diffraction  and  x-ray  photoelectron 
spectroscopy.  The  films  showed  no  crystalline  structure  in  our  x-ray  diffraction 
measurements.  X-ray  photoelectron  spectroscopy  of  the  samples  showed  no 

discernable  shift  in  the  Ti-O  binding  energy  between  slightly  absorbing  and 
nonabaorbing  films.  This  indicates  that  the  changes  affecting  optical  absorption  are 
small. 

From  the  above  data  we  propose  the  following  model:  As  we  increase  the 
number  of  arriving  oxygen  ions,  the  number  of  oxygen  vacancies  decreases.  As 
mentioned  earlier,  TiOj  is  known  to  be  oxygen  deficient  when  thermally  evaporated. 
These  oxygen  vacancies  are  filled  by  chemically  active  oxygen  ions  absorbed  at  the 
substrate.  This  simultaneously  increases  the  refractive  index  and  decreases  the 
extinction  coefficient. 

The  filling  of  oxygen  vacancies  is  a  balance  between  opposing  factors.  This  is 
demonstrated  by  examining  the  optical  constants  as  a  function  of  7.  As  7  increases 

from  0,  we  fill  oxygen  vacancies  until  with  0.40  £  7  £  0.45.  the  vacancies  have 
been  saturated  and  the  preferential  sputtering  of  oxygen  out  of  the  film  becomes 
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significant.  Although  the  sputtering  yield  is  low  for  300  eV  ions  at  7  £  0.5,  the 
statistics  begin  to  favor  the  preferential  sputtering  of  oxygen  from  the  surface  of  the 
growing  film  because  the  peak  ion  flux  is  so  large  and  the  Ti-0  chemical  bond  in 


TiO,  is  weak.11 


Sputtering,  a  function  of  ion  energy,  should  increase  at  higher  ion 
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energy.  Indeed,  we  noted  that  our  films  were  damaged  at  a  lower  7  for  500  eV 
intermittent  ion  bombardment.  This  apparently  occurs  when  the  population  of  active 
oxygen  at  the  substrate  is  below  that  needed  to  satisfy  the  oxygen  deficiency. 

In  conclusion,  intermittent  km  bombardment  can  increase  the  refractive  index 
and  reduce  the  absorption  of  TiO,  films.  The  range  of  refractive  indices  reported 
here  is  similar  to  that  reported  elsewhere9  for  LAD  without  substrate  rotation.  This 
demonstrates  that  IAD  can  be  considered  for  production  environments. 

The  authors  would  like  to  acknowledge  the  support  of  the  Air  Force  Office  of 
Scientific  Research  through  the  University  Research  Initiative  Program. 
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Figure  Captions 

Figure  1.  Mean  refractive  index  at  550  run  of  TiOj  films  bombarded  with  300  eV 
oxygen  ions  as  a  function  of  the  average  number  of  ions  per  arriving 
molecule,  o  -  38  cm  from  ion  gun  to  optic.  A  -  30  cm  from  ion  gun 
to  optic. 

Figure  2.  Extinction  coefficient  at  550  nm  of  TiO,  films  bombarded  with  300  eV 
oxygen  ions  as  a  function  of  the  average  number  of  ions  per  arriving 
molecule,  o  -  38  cm  from  ion  gun  to  optic.  A  -  30  cm  from  ion  gun 
to  optic. 

Figure  3.  Relative  inhomogeneity  at  550  nm  of  TiO,  films  bombarded  with  300  eV 
oxygen  ions  as  a  function  of  the  average  number  of  ions  per  arriving 
molecule,  o  -  38  cm  from  ion  gun  to  optic.  A  -  30  cm  from  ion  gun 
to  optic. 

Figure  4.  Vacuum-to-air  shift  at  450  nm  of  TiO,  films  bombarded  with  300  eV 
oxygen  ions  as  a  function  of  the  average  number  of  ions  per  arriving 
molecule,  o  -  38  cm  from  ion  gun  to  optic.  A  -  30  cm  from  ion  gun 
to  optic. 
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ABSTRACT 


Aluminum  nitride  thin  films  have  been  deposited  by  ion-assisted  deposition.  Aluminum 
was  electron-beam  evaporated  onto  substrates  with  simultaneous  nitrogen  ion  bombardment 
Rutherford  backscattering  spectrometry  showed  that  nitrogen-to-aluminum  ratios  of  one  or 
greater  could  be  achieved  with  sufficient  nitrogen  ion  fluxes.  This  excess  nitrogen  apparently 
degrades  the  optical  properties  of  the  films  in  the  visible.  Annealing  at  500°C  improves  the 
optical  properties  drastically  at  the  expense  of  a  slight  oxygen  diffusion  into  the  films.  Finally, 
aluminum  oxynitride  films  were  deposited  by  adding  an  oxygen  backfill  to  the  vacuum 
chamber  during  deposition.  These  films  had  very  similar  optical  properties  to  the  annealed 
nitride  films. 


INTRODUCTION 


Aluminum  nitride  (AIN)  is  a  very  attractive  coating  material  for  many  applications.  Its 
mechanical  strength  and  chemical  stability  have  resulted  in  its  use  aw  a  protective  layer  on 
machine  tools  and  as  an  insulating  layer  in  microelectronic  applications.  More  recently,  its 
relatively  wide  bandgap  (6.5  eV)  has  led  to  its  consideration  as  an  optical  coating  material  for 
use  in  the  visible  and  near-ultraviolet  regions  (1,2). 

Aluminum  nitride  films  have  been  deposited  by  many  methods,  including  chemical-vapor 
deposition  (31,  reactive  ion  sputtering  (41  ion  implantation  [5],  and  dual-beam  ion  sputtering  [6,71 
In  this  work  we  have  bomarded  a  growing  aluminum  film  with  a  nitrogen  ion  beam.  With  the 
aluminum  thermally  evaporated  from  an  electron-beam  source.  It  was  hoped  that  the  use  of  a 
thermal  source  could  provide  a  large  aluminum  deposition  rate  and  allow  the  deposition  of 
AIN  over  a  wide  area  in  a  system  where  only  one  ion  source  was  available.  Martin  and 
Netterfield  have  recently  reported  the  successful  ion-assisted  deposition  (IAD)  of  Si3N4  [81 
suggesting  that  aluminum  nitride  deposition  should  also  be  feasible  by  IAD.  We  have  also 
deposited  aluminum  oxynitride  films  by  adding  oxygen  to  the  chamber  backfill  during  the 
depositions. 


DEPOSITION  CONDITIONS 


The  films  were  deposited  in  a  Balzers  BAK  760  box  coater  using  the  geometry  shown  in 
Figure  1.  The  substrates,  glass,  fused  silica,  and  graphite,  were  selected  for  their  compatibility 
with  various  forms  of  analysis  and  were  held  stationary  during  the  deposition.  All  coating  runs 
were  performed  at  100°C  with  a  residual  pressure  of  approximately  6x1 0‘7  torr.  Aluminum 
starting  material  of  99.999%  purity  was  evaporated  by  an  electron  beam  source  from  an 
intermetallic  crucible  at  an  aluminum  deposition  rate  of  0.2  to  0.7  nm/s.  The  aluminum 
deposition  rate  was  measured  by  a  quartz  crystal  rate  monitor  shielded  from  the  ion  beam.  In 
most  of  the  depositions,  the  chamber  was  filled  with  6x1  O'3  torr  of  N2  gas  through  the  ion 
gun.  Gas  could  also  be  bled  into  the  chamber  through  a  leak  valve  in  the  side  of  the 
chamber,  bypassing  the  plasma  discharge  in  the  ion  gun.  The  substrates  were  ion  precleaned 
before  deposition  and  bombarded  during  deposition  by  nitrogen  ions  from  a  3-cm-aperture 
Kaufman  ion  source  equipped  with  tungsten  filaments  and  graphite  grids.  The  filament 
neutralizer  emitted  thermal  electrons  to  neutralize  the  beam,  preventing  charge  buildup  on  the 
dielectric  substrates.  A  Faraday  cup  measured  the  ion  current  density.  Aluminum  nitride  films 
were  deposited  by  bombarding  the  growing  aluminum  films  with  a  nitrogen  ion  beam  with  ion 
energies  of  250,  500,  750,  or  1000  eV.  51 
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Figure  1.  Schematic  of  Balzers  BAK  760  coating  plant  configured  for  ion- 
assisted  deposition. 


RUTHERFORD  BACKSCATTERING  SPECTROMETRY 


Film  stoichiometries  can  be  measured  very  accurately  (to  better  than  one  percent)  using 
Rutherford  backscattering  spectrometry  (RBS)  [9J.  Graphite  substrates  were  used  in  this  work 
because  the  silicon  and  oxygen  RBS  signals  from  glass  substrates  overlap  those  from 
aluminum,  nitrogen,  and  oxygen  in  the  films,  decreasing  measurement  accuracy  and 
complicating  the  analysis.  In  addition,  aluminum  nitride  films  must  be  less  than  100  nm  thick  to 
prevent  the  oxygen  and  nitrogen  backscattering  peaks  from  overlapping  for  1.892  MeV  He* 
ions. 

We  first  deposited  an  aluminum  film  at  a  rate  of  02  nm/s  in  a  nitrogen  backfill  of 
6x10"5  torr  to  determine  if  (a)  this  nitrogen  would  be  adsorbed  into  the  films  and  (b)  whether 
our  base  vacuum  was  sufficiently  free  of  oxygen  and  water  vapor  to  prevent  significant 
oxidation  of  the  films.  The  nitrogen  content  of  these  films  was  less  than  2  at  %.  The  bulk  of 
the  film  contained  1.1  ±  02  at  %  oxygen,  with  greater  oxidation  observed  near  both  the 
substrate-film  and  film-air  interfaces. 

Table  1  shows  the  stoichiometries  of  aluminum  nitride  films  as  determined  by  RBS.  The 
stoichiometries  are  expressed  in  the  form  AINxOy.  This  is  simply  a  convenient  notation,  and 
does  not  imply  the  formation  of  a  compound  with  the  indicated  stoichiometry.  The  ratio  7  of 
incident  nitrogen  atoms  to  aluminum  atoms  at  the  substrate  has  been  calculated  under  the 
assumption  that  the  ion  beam  consists  of  N2*  ions. 

Most  of  the  films  have  a  stoichiometry  between  AINqs  and  AIN  1  2.  with  a  large  number  of 
the  films  having  superstoichiometric  nitrogen  concentrations  (Figure  2).  Harper,  Cuomo,  and 
Hentzeil  have  previously  reported  [6]  that  dual-beam  sputtered  films  with  a  similar  energy  for 
the  ion  beam  bombarding  the  substrate  are  stoichiometric  if  the  nitrogen-to-aluminum  arrival 
rate  at  the  substrate  is  greater  than  one.  They  hypothesize  that  this  self-limiting  effect  in  the 
stoichiometry  is  due  to  diffusion  of  nitrogen  to  the  surface  of  the  films  followed  by  desorption 
of  N2  from  the  surface.  We  observed  no  such  self-limiting  effect  with  ion-assisted  deposition. 
A  1:1  ratio  is  achieved  only  for  particular  sets  of  deposition  conditions  rather  than  for  whole 
regions  of  the  deposition  parameter  space  as  observed  with  dual-beam  sputtering.  Since  the 
only  significant  difference  between  the  two  methods  is  the  aluminum  source,  it  is  possible  that 
highly  energetic  sputtered  aluminum  atoms  are  responsible  for  the  diffusion  of  excess  nitrogen 
atoms  to  the  surface  of  the  films  where  they  could  be  desorbed.  Bombardment  by  the 
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Table  1.  Stoichiometries  of  AIN  films  (AINxOy). 


A!  dep.  rate 
(nm/sec) 


Ion  Beam 

Energy  Current  Density 
(eV)  (/iA/cm^) 


90  £  20 
100 
140 
130 
190 


Gamma 

(N  atoms/AI  atom) 


0.9  t  0.3 
1.0  t  0.4 
1.5  t  0.5 
1.1  £  0.4 
2.0  £  0.5 

0.5  £  0.1 
0.8  £  02 
0.9  t  02 

1.0  £  0.3 
2.0  t  0.5 
22  t  0.5 
02  £  02 
0.4  £  0.1 
0.9  t  02 
1.3  £  02 
0.9  t  02 

1.8  £  0.5 
2.7  £  0.6 


Stoichiometry 


0.89  *  0.01 
1.01  £  0.01 
121  £  0.02 
0.81  £  0.02 
124  £  0.02 

0.58  £  0.01 
1.04  £  0.02 
1.15  £  0.02 

1.03  £  0.02 
123  £  0.03 
1.17  £  0.02 
0.82  £  0.01 
0.55  £  0.01 
127  £  0.02 
1.11  £  0.02 
126  £  0.02 

129  £  0.03 
1.58  £  0.10 


0.04  £  0.006 
0.03  £  0.01 
0.04  £  0.01 
0.06  £  0.01 
0.03  £  0.003 

0.03  £  0.01 
0.03  £  0.01 
0.024  £  0.006 

0.05  £  0.01 
0.12  £  0.03 
0.09  £  0.02 
0.02  £  0.01 
0.023  £  0.004 
0.03  £  0.003 
0.08  £  0.02 
0.027  £  0.003 

021  £  0.01 
0.52  £  0.04 


nitrogen  beam  alone  presumably  does  not  provide  sufficient  mobility  to  the  growing  film  to 
allow  sufficient  nitrogen  diffusion. 

The  RBS  results  also  show  that  the  films  have  very  high  packing  densities:  The  majority 
of  the  films  show  no  measurable  amount  of  oxygen  below  the  surface,  suggesting  that  any 
oxygen  adsorbed  from  the  residual  chamber  atmosphere  during  deposition  is  resputtered,  and 
that  the  films  are  too  dense  for  water  to  migrate  through  the  films  after  venting  the  system  to 
atmosphere.  There  is  an  oxidized  layer  at  the  surface  of  the  films  with  a  thickness  of 
15  £  7.5  nm.  The  oxygen-to-aluminum  ratios  y  in  Table  1  are  averages  across  the  film, 
including  this  surface  layer. 

A  final  observation  based  on  the  RBS  analysis  is  that  the  film  deposition  rate  can  be 
increased  if  a  sufficient  ion  flux  is  available.  We  have  shown  that  300  m A/cm^  of  750  eV  ions 
is  sufficient  to  produce  a  superstoichiometric  film  for  an  aluminum  deposition  rate  of 
0.65  nm/s. 

The  RBS  results  indicate  several  sets  of  of  deposition  conditions  which  yield  a  1:1 
aluminum-to-nitrogen  ratio.  For  example,  evaporating  aluminum  at  a  rate  of  0.2  nm/s  while 
bombarding  with  a  750-eV,  lOO-^A/cm^  ion  beam  should  result  in  a  stoichiometric  film. 
However,  we  shall  show  in  the  next  section  that  films  deposited  with  parameters  similar  to  this 
are  optically  absorbing 

OPTICAL  MEASUREMENTS 

Transmittance  and  reflectance  measurements  have  been  performed  on  films  deposited 
onto  fused  silica  substrates  in  a  Cary  14  dual-beam  spectrophotometer  over  a  wavelength 
range  of  02  to  2.0  fim.  Films  deposited  with  an  insufficient  ion-beam  flux  are  opaque  with  a 
shiny  or  dull  grey  appearance  depending  on  the  deposition  conditions.  The  sheet  resistance 
of  these  films  has  been  measured  with  a  four-point  probe,  which  shows  that  the  resistance  of 
films  rises  above  the  3000  Q/square  upper  limit  of  our  instrument  when  the  films  begin  to 
show  any  transmittance  in  the  visible.  As  the  film  approach  the  correct  stoichiometry,  the 
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transmittance  increases.  In  even  the  best  case,  however,  the  films  begin  to  absorb  at 
wavelengths  shorter  than  500  nm  (Rgure  3).  Holmes  has  observed  a  similar  behavior  with 
reactive  ion  beam  sputtered  films  [4],  which  he  attributed  to  nitrogen  vacancies  or  defect 
states,  although  his  films  showed  no  absorption  until  450  nm.  Our  films  are  in  general 
nonabsorbing  at  wavelengths  longer  than  1  pm,  and  we  have  therefore  used  the  quarter-wave 
transmittance  in  the  1  to  2  Mm  region  to  calculate  a  refractive  index  for  these  films.  As  the 
dispersion  curve  of  A}N  is  essentially  flat  over  this  wavelength  region,  calculations  performed 
at  slightly  different  wavelengths  are  still  comparable.  Table  2  gives  the  refractive  indices  of 
these  films;  these  values  are  slightly  lower  than  the  literature  values  of  1.90  to  1.95  [1]. 


0  t  2  3  WAVELENGTH  (••) 

V 


Figure  2.  Stoichiometry  of  aluminum  nitride  Figure  3.  Transmittance  of  typical  AIN 

films,  expressed  in  the  form  AINX,  as  a  films  deposited  by  ion-assisted  deposition, 

function  of  the  ratio  7  of  nitrogen  to 
aluminum  arrival  rates. 


Table  2.  Near  infrared  refractive  indices  of  AIN  films. 


Ion  Beam 


Al  dap.  rate 
(nm/s) 

Energy 

(eV) 

Current  Density 
(jj.A/cm£) 

Gamma 

(N  atoms/AI  atom) 

Refractive 
index  n 

02 

750 

185  ±  20 

1.9  ±  0.6 

1.87  ±  0.02 

0.3 

500 

150 

1.0  t  0.3 

1.89 

0.3 

500 

190 

1.3  t  0.3 

1.82 

0.4 

750 

290 

1.5  ±  02 

1.86 

0.6 

750 

400 

1.4  i  0.2 

1.84 

ANNEAUNG  STUDIES 

In  an  attempt  to  improve  the  optical  properties  of  the  films,  a  number  of  films  were 
annealed  in  a  tube  furnace  at  500°C.  Dry  nitrogen  was  flowed  through  the  tube  to  reduce 
the  oxygen  and  water  content  in  the  furnace,  but  the  tube  was  not  sealed.  Table  3  shows  the 
effect  of  annealing  on  three  RBS  samples  with  widely  varying  stoichiometries.  We  see  that 
the  nitrogen  content  of  the  films  remains  the  same  or  decreases  after  210  minutes,  and  that  the 
oxygen  content  of  the  films  increases  slightly.  The  oxygen  in  the  original  films  was 
concentrated  in  a  surface  layer  which  became  thicker  upon  annealing  as  oxygen  penetrated 
into  the  films. 

The  annealing  of  thicker  samples  on  fused  silica  substrates  showed  that  the  optical 
properties  of  the  films  were  improved  substantially  (Figure  4),  with  no  absorption  observed  at 
wavelengths  longer  than  500  nm,  and  the  absorption  edge  moved  below  350  nm.  Additional 


Table  3.  Effect  of  annealing  on  stoichiometry  of  AJNxOy  films. 


ion  Beam 


Al  dep. 
rate 

(nm/s) 

Energy 

(eV) 

Current 

Density 

0*A/cm2) 

0.45 

750 

100 

025 

250 

130 

0.65 

750 

290 

Stoichiometry 


Unanneaied 


0.55  *  0.01 

0.81  t  0.02 
1.26  t  0.02 


0.023  t  0.004 
0.06  ±  0.01 
0.027  ±  0.003 


Annealed 


0.56  ±  0.01 
0.74  *  0.02 
1.15  ±  0.02 


0.068  t  0.004 
0.15  t  0.01 
0.054  t  0.005 


annealing  at  500°C  did  not  affect  the  transmittance  of  the  films,  although  additional  oxygen  is 
certainly  diffusing  into  the  films.  This  suggests  that  unsatisfied  bonds  in  the  films  are  being 
filled  quickly  by  oxygen  atoms  or  complexes,  eliminating  the  dangling-bond  related  absorption. 
The  cause  of  the  remaining  absorption  is  still  not  positively  identified. 

OXYNITRIDE  DEPOSITION 

Mixed  compounds  of  oxides  and  nitrides  are  useful  in  tailoring  the  refractive  index  of 
coating  materials.  Aluminum  oxynitrides  should  span  the  refractive  index  range  between  1.65 
and  1.9  in  the  visible.  We  have  deposited  aluminum  oxynitride  films  by  adding  3x10*5  torr  of 
oxygen  to  the  vacuum  chamber  during  the  deposition.  This  oxygen  then  occupies  adsorption 
sites  created  by  the  nitrogen  ion  bombardment  Figure  5  shows  the  transmittance  of  two  films 
deposited  with  an  aluminum  deposition  rate  of  0.6  nm/sec  and  bombardment  by  750  eV  N2* 
ions.  Both  films  have  significantly  better  optical  properties  than  the  aluminum  nitride  films 
discussed  above.  This  is  presumably  due  to  adsorbed  oxygen  atoms  filling  vacancies  created 
by  the  ion  bombardment  The  films  still  show  an  absorption  edge  more  typical  of  the  nitride 
materiat  The  200  pA/crc?-  sample  has  a  refractive  index  of  1.93  while  the  300  pA/cm^  film 
has  an  index  of  1.82  because  of  its  greater  oxygen  content  The  refractive  index  can  thus  be 
varied  by  changing  the  ion  beam  current  density.  An  alternative  method  would  be  to  pass 
oxygen  through  the  ion  gun.  This  oxygen  would  be  more  effectively  implanted  into  the  films, 
resulting  in  a  larger  oxygen  content  for  the  same  ion  beam  current  density.  The  refractive 
index  of  1.93  is  surprisingly  high,  and  may  suggest  that  the  amount  of  adsorbed  oxygen  in  the 
film  is  quite  small.  1 00  1 — - - - - - - 
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Figure  4.  Transmittance  of  an  AIN  film  Figure  5.  Transmittance  of  aluminum 

deposited  by  I  AD  before  and  after  annealing  oxynitride  films  deposited  by  l  AD. 

for  210  minutes  at  50CPC. 


CONCLUSIONS 


Aluminum  nitride  thin  films  can  be  deposited  by  ion-assisted  deposition;  however,  the 
diffusion  and  desorption  of  excess  nitrogen  observed  for  dual-beam  ion  sputtering  is  not 
evident  with  ion-assisted  deposition.  An  excess  of  nitrogen  is  therefore  present  in  many  of  the 
films.  A  nitrogen  deficiency  results  in  opaque  conducting  films,  while  a  nitrogen  excess  results 
in  insulating  films  with  absorption  in  the  visible  region  of  the  spectrum.  This  absorption  can 
be  decreased  by  annealing  in  a  predominantly  nitrogen  atmosphere  at  500°C,  suggesting  that 
the  absorption  is  due  to  nitrogen  vacancies.  Absorption  is  still  present  below  500  nm,  but  the 
films  now  possess  a  respectable  transmittance  in  the  near  ultraviolet 

Adding  an  oxygen  backfill  to  the  chamber  during  deposition  allows  the  deposition  of 
aluminum  oxynitride  films.  These  films  have  significantly  better  optical  properties  than  the 
nitride  films  and  allow  the  deposition  of  films  with  a  refractive  index  between  those  of  AI2O3 
and  AIN. 
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Arjomandy,  M.  F.  Hinedi,  J.  Oder,  J.  Serveld,  and  G.  van  Zyll  for  assistance  with  RBS  data 
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STOICHIOMETRY  OF  HIGH  Tc  THIN  FILMS 
BY  HIGH  ENERGY  BACKSCATTERING 
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TUCSON,  AZ  85721 


ABSTRACT 


ABSOLUTE  AREAL  DENSITIES 
(ATOMS/SQ  CM)  AND  RELATIVE 
AREAL  DENSITIES  OF  HIGH  Tc  THIN 
FILMS  MAY  BE  DETERMINED  BY  HIGH 
ENERGY  BACKSCATTERING,  USING  3 
TO  5  MEV  He  ANALYSIS  BEAMS,  TO 
ACCURACIES  OF  ABOUT  3%  AND 
BETTER  THAN  \%,  RESPECTIVELY. 

WE  SHOW  THREE  EXAMPLES  OF 
THE  ANALYSIS  OF  YBaCuO  FILMS  IN 
THE  THICKNESS  RANGE  2500A  TO 
9000A. 


FIGURE  1 

ORDINARY  1 .9  MEV  He  RBS  SPECTRUM 
OF  A  THIN  YBaCuO  FILM  (SAMPLE  A) 
HEAVILY  CONTAMINATED  WITH  Al 
AND  W.  NOTE  THE  INTERFERENCE  OF 
THE  Cu,  Y,  Ba  AND  W  SIGNALS. 
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FIGURE  2 


HIGH  ENERGY  3.8  MEV  BACK- 
SCATTERING  SPECTRUM  OF  THE 
SAMPLE  A.  ONLY  THE  Ba  AND  W 
SIGNALS  INTERFERE  IN  THIS 
SPECTRUM.  THE  He-0  SCATTERING 
IS  NON-RUTHERFORD. 
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FIGURE  3 


COMPUTER  DECONVOLUTION  OF  THE 
INTERFERING  Ba  AND  W  SIGNALS 
SHOWN  IN  FIGURE  2. 
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TABLE  1 


RESULTS  FOR  SAMPLE  A 
ELEMENT  Nt  ATOMIC  X 


Y 

0.51+0.02 

2.4+0. 1 

Ba 

0.33+0.02 

1. 5+0.1 

Cu 

1.67+0.05 

7.7+0. 1 

0 

10.0+0.4 

46.0+0.6 

Al 

9.03+0.29 

41.7+0.6 

W 

0.083+0.005 

0.4+0. 1 

Ar 

0.069+0.010 

0.3+0. 1 

HI  =  AREAL  DENSITY  (El 7  ATOMS/SQ  CM) 
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LIMITS  ON  THE  ACCURACY  OP  STOICHIOMETRY  DETERMINED  BY 
RUTHERFORD  BACKSCATTERING  USING  COMPUTER  PEAK  FITTING. 


L.  C.  MCINTYRE.  JR..  M.  D.  ASHBAUGH  AND  J.  A.  LEAVITT 
Department  of  Physics,  University  of  Arizona.  Tucson.  AZ  8S721 


ABSTRACT 


Rutherford  Backscattering  Spectrometry  (RBS)  using  MeV  ions  is  capable  of  routinely 
measuring  absolute  atomic  areal  densities  to  an  accuracy  of  about  3%  and  relative  areal 
densities  to  a  fraction  of  1%  if  the  elastic  scattering  peaks  from  different  elements  do  not 
overlap.  The  accuracy  of  areal  density  measurements  can  be  seriously  reduced  in  cases 
where  the  sample  thickness  is  large  enough  to  cause  overlapping  peaks.  We  report  here  an 
investigation  of  the  use  of  standard  computer  fitting  techniques  to  analyze  overlapping  peaks 
in  elastic  scattering  spectra  and  its  effect  on  the  accuracy  of  the  final  results  of  an  RBS 
analysis. 


INTRODUCTION 


Rutherford  Backscattering  Spectroscopy  using  MeV  energy  beams  of  He  and  other 
ions  has  proved  to  be  a  useful  and  accurate  method  for  measuring  elemental  areal  densities 
and  stoichiometry  of  thin  films  and.  in  certain  cases,  of  surface  layers  of  bulk  materials. 
The  absolute  accuracy  of  this  method  is  limited  to  about  3%.  due  mainly  to  uncertainties  in 
charge  collection  of  the  incident  beam  and  the  detector  solid  angle.  These  quantities  are 
normally  calibrated  with  a  standard  Bi  in  Si  target  [1]  in  which  the  Bi  areal  density  is 
known  to  at  least  this  accuracy.  Stoichiometry,  or  relative  atomic  concentrations,  can  be 
determined  more  accurately  since  the  factors  mentioned  above,  which  contribute  the  bulk  of 
the  uncertainty,  cancel  in  these  calculations. 

Difficulties  arise  in  cases  where  the  films  contain  similiar  mass  elements  and/or  are  of 
sufficient  thickness  to  cause  overlapping  of  the  elemental  (or  in  some  cases  isotopic)  peaks  in 
the  backscattering  energy  spectrum.  Several  computer  programs  have  been  reported  [2-6] 
that  fit  complex  backscattering  spectra  and  determine  areal  densities  by  modeling  the 
spectrum  using  the  parameters  of  the  experiment,  including  detector  resolution,  and  the 
stopping  cross  section  of  the  incident  ion  in  the  film.  We  present  here  a  somewhat  different 
and  simpler  approach  which  assumes  each  element  in  a  single  layer  film  contributes  a 
standard  shape  to  the  spectrum,  and  the  complex  spectrum  is  a  sum  of  several  of  these 
standard  shapes.  The  aim  is  simply  to  deconvolve  several  of  these  overlapping  shapes  and 
determine  the  individual  areas.  No  experimental  parameters  or  stopping-power  information 
is  used  in  the  program.  The  only  input  information  required  is  an  initial  estimate  of  the 
parameters  of  the  standard  shapes.  The  output  is  the  best-fit  parameters  from  which  peak 
areas  can  be  calculated.  We  have  found  this  procedure  to  be  useful  and  efficient  in  routine 
RBS  anaylsis  using  a  standard  IBM  PC  microcomputer. 


FITTING  FUNCTION  AND  PROCEDURE 


A  thin  film  consisting  of  a  uniform  distribution  of  a  single  isotope,  which  is  either 
self-supporting  or  deposited  on  a  substrate  containing  atoms  less  massive  than  those  in  the 
film,  produces  a  backscattering  energy  spectrum  with  a  shape  similiar  to  that  shown  in 
Fig.  1. 

We  have  written  a  piecewise  analytic  function  with  six  parameters  to  represent  this 
shape.  It  is  a  trapezoid  with  a  half  Gaussian  on  each  end.  The  parameters  are  shown  in 
Fig.  1  and  represent  the  position  and  height  of  each  side  of  the  trapezoid  and  the  width  of 
each  Gaussian.  The  procedure  used  to  deconvolve  several  of  these  overlapping  shapes 
caused  by  multiple  elements  in  an  actual  spectrum  is  to  use  a  nonlinear  least-squares  search 
routine  to  find  the  shape  parameters  (six  per  element)  which  best  reproduce  the  data. 
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Fig.  1.  Piecewise  anaiyic  shape  representing  the  backscattering  spectrum 
from  a  single  element.  The  numbers  1-6  indicate  the  parameters  of  the  shape. 

We  have  used  the  subroutine  CURFIT  in  Ref.  [7]  to  determine  these  parameters.  The 
program  was  written  in  IBM  Professional  FORTRAN  and  is  run  on  a  standard  IBM  PC 
microcomputer  (with  an  8087  math  coprocessor).  Some  early  results  using  this  procedure 
have  been  reported  previously  [8]. 

This  method  was  first  tested  on  a  fully  resolved,  single-element  peak  to  assure  that 
our  analytic  function  would  adequately  represent  the  data.  When  this  proved  satisfactory, 
we  proceeded  to  examine  two-  and  three-element  cases  with  various  degrees  of  overlap. 
Some  of  these  cases,  together  with  comparisons  to  data  taken  at  higher  energy  where  the 
elemental  peaks  are  resolved,  will  be  discussed  in  the  next  section.  We  have  applied  this 
method  to  a  number  of  special  cases  in  which  we  invoke  particular  relationships  between 
the  parameters  with  some  degree  of  success.  If  any  of  the  constituents  has  a  nonuniform 
depth  distribution,  this  method  is  not  appropriate.  We  have  not  corrected  for  the  existence 
of  several  isotopes  for  a  given  element;  however,  this  refinement  could  be  added. 

TEST  CASES 

Tantalum  Molybdenum  Superlattice  Film 

An  example  of  the  fitting  procedure  described  above  is  shown  in  Fig.  2.  which  shows 
the  backscattering  spectrum  of  1892  keV  He  ions  scattered  from  a  TaMo  superlattice  film 
approximately  300  nm  thick  on  a  carbon  substrate. 

At  this  incident  energy,  the  peaks  from  the  two  elements  are  not  resolved.  The  solid 
curve  in  Fig.  2  is  the  result  of  a  10-parameter  fit  to  the  data  using  two  of  our  standard 
shapes.  We  used  10  instead  of  12  parameters  by  setting  the  two  left  Gaussian  widths  and 
the  two  right  Gaussian  widths  equal.  The  final  parameters  were  obtained  with  five 
iterations  of  the  CURFIT  procedure,  which  took  about  30  sec  per  iteration. 

Separate  peak  areas  and  their  uncertainties  (based  on  the  uncertainties  in  the 
parameters  given  by  CURFIT)  were  calculated  from  the  final  parameters,  and  the 
stoichiometry  and  areal  density  of  each  element  was  calculated.  This  same  film  was  then 
analyzed  using  4700  keV  He  ions,  which  resulted  in  a  backscattering  spectrum  in  which  the 
two  peaks  are  completely  resolved.  Areal  densities  and  stoichiometry  were  calculated  from 
the  resolved  peaks;  the  two  sets  of  results  are  shown  in  Table  I. 

We  see  that  the  areal  densities  measured  in  the  two  cases  are  statistically  consistent  and 
have  the  same  uncertainty.  This  implies  that  the  3%  uncertainty  in  charge  collection  and 
solid  angle  dominates  any  additional  uncertainties  produced  by  the  computer  deconvolution 
in  the  low-energy  case.  The  uncertainty  in  atomic  fractions,  however,  is  a  factor  of  five 
worse  in  the  computer  fitted  case. 
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Fig.  2.  Backscattering  spectrum  of  1892  keV  He  ions  on  a  TaMo  superlattice 
approximately  300  am  thick.  The  solid  curve  is  a  10- parameter  fit  to  the  data. 


Table  I.  Areal  density  and  atomic  fractions  of  TaMo  film  (Nt  is  atoms/cm2  in  units  of  101T). 


1892  keV  Computer  Fit 

4700  keV  Resolved  Peaks 

Nt 

F 

Nt 

F 

Mo 

9.2  ±  0.3 

0.497  t  0.005 

9.3  t  0.3 

0.501  t  0.001 

Ta 

9.3  ±  0.3 

0.503  t  0.005 

9.2  ±  0.3 

0.499  ±  0.001 

Aluminimum  Oxide-< Argon)  Film 

An  example  with  three  overlapping  peaks  is  shown  in  Fig.  3.  This  is  a  spectrum 
obtained  by  backscattering  1892  keV  He  ions  from  an  aluminimum  oxide  film  with  nominal 
thickness  about  600  nm  containing  argon.  The  solid  curve  is  a  14-parameter  fit  to  the  data 
using  three  of  our  standard  shapes.  As  in  the  two-peak  case,  the  three  left  as  well  as  the 
three  right  Gaussian  widths  were  constrained  to  be  equal. 

This  film  was  also  analyzed  using  3776  keV  He  ions,  which  produced  a  spectrum  in 
which  the  three  peaks  were  completely  resolved.  Table  II  shows  a  comparison  of  the  results 
o burned  in  the  two  cases. 

We  see  again  that  the  two  results  are  consistent,  the  areal  densities  have  similiar 
uncertainties,  and  the  uncertainties  in  the  atomic  fractions  are  considerably  better  in  the 
separated  peak  case.  It  should  be  noted  that  a  correction  [9]  was  applied  to  the  oxygen 
peak  in  the  3776  keV  case  to  account  for  non-Rutherford  scattering. 


SUMMARY  AND  CONCLUSION 

We  have  described  a  simple  and  rapid  application  of  standard  methods  for 
deconvolving  overlapping  peaks  in  an  elastic-scattering  spectrum  using  a  microcomputer. 
Its  use  is  limited  to  films  with  uniformly  distributed  elemenul  constituents  within  a  given 
layer.  We  have  found  this  procedure  to  be  convenient  and  useful  in  routine  RBS  analysis, 
and  it  requires  no  experimental  parameters  or  ion-stopping  cross  sections.  The  resulting 
individual  peak  areas  are  subsequently  included  in  the  total  analysis  of  the  film  together 
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Fig.  3.  Backscattering  spectrum  of  1892  keV  H«  ions  oa  a  AJ,0,  film 
approximately  600  am  thick  coatainiag  Ar.  The  solid  curve  is  a 
14-parameter  fit  to  the  data. 


Table  II.  Areal  density  and  atomic  fractions  of  Al,Oj-Ar  film  (Nt  is  atoms/cm2 
in  units  of  101*). 


1892  keV  Computer  Fit  3776  keV  Resolved  Peaks 


Nt 

F 

Nt 

F 

A1 

2.65  ±  0:09 

0.368  ±  0.005 

2.63  t  0.08 

0.367  ±  0.003 

0 

4.16  ±0.16 

0.577  ±  0.006 

4.12  ±0.13 

0.575  ±  0.003 

Ar 

0.40  t  0.02 

0.055  ±  0.001 

0.41  s  0.01 

0.058  ±  0.001 

with  any  resolved  peak  areas.  The  examples  given  here  involve  thick  targets  with  resulting 
wide  peak  shapes;  however,  in  the  majority  of  cases  we  have  encountered  involving  thinner 
films,  a  typical  iteration  time  is  5-15  sec.  and  the  average  number  of  iterations  required  is 
about  five.  We  see  that  areal  densities  with  comparable  accuracy  can  be  obtained  from 
spectra  containing  resolved  or  overlapping  peaks,  but  the  accuracy  of  the  resulting 
stoichiometry  is  still  considerably  greater  in  the  resolved  peak  cases. 

We  would  like  to  thank  C.  Falco.  J.  Makous,  and  L.  Gignal  for  providing  us  with  the 
films  used  in  the  examples  presented  here.  We  also  thank  P.  Stoss.  B.  Dezfouly-Arjomandy. 
M.  F.  Hinedi.  J.  Oder.  J.  Seerveld.  and  G.  van  Zyll  for  help  in  taking  the  data.  This  work 
was  partially  supported  by  IBM-Tucson.  The  Optical  Data  Storage  Center,  and  the  Air  Force 
Office  of  Scientific  Research  through  the  University  Research  Initiative  Program. 
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ELASTIC  PROPERTIES  OF  EPITAXIAL  ZnSe  (001)  FILMS  ON  GaAs 
MEASURED  BY  BRILLOUIN  SPECTROSCOPY 
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Brillouin  scattering  experiments  have  been  performed  on  a  series  of  epitaxially 
grown  ZnSe(00l)  films  on  GaAs  (001)  of  thicknesses  between  0.23  fim  and  2.0  fim.  No 
thickness-dependent  modifications  of  the  elastic  constants  have  been  found.  The 
obtained  room  temperature  values  of  Cu  -  87.0,  Cl2  -  54.7  and  Cu  -  39.1  GPa  are  in 
good  agreement  with  literature  values  of  ZnSe  bulk  material.  The  calculation  of  the 
theoretical  cross  sections  shows  a  good  agreement  with  the  experimental  data. 


PACS  numbers:  78.3S.+C,  62.20. Dc 

^  Also  at  Gruppo  Nazionale  Struttura  della  Materia,  Universita  di  Modena,  Italy 


Epitaxial  films  of  ZnSe  have  drawn  considerable  interest  for  their  physical 
properties,  as  well  as  their  applications  in  optical  and  ultrasonic  devices.  Since  the 
ZnSe  lattice  parameter  exceeds  that  of  GaAs  at  room  temperature  by  0.27%,  epitaxially 
grown  ZnSe(001)  films  mi  GaAs(001)  are  compressively  strained  in  the  plane  of  the 
interface,  and  because  of  the  resulting  tetragonal  distortion  of  the  lattice  (Poisson  effect), 
they  are  tensionally  strained  normal  to  this  plane.  Mitsuhashi  et  al.1  have  measured  by 
X-ray  diffraction  the  dependence  of  the  perpendicular  lattice  parameter  of  ZnSe  on  the 
layer  thickness  for  films  grown  on  GaAs  (001)  by  metal-organic  chemical  vapor 
deposition  (MOCVD)  at  400  °C.  They  find  that  layers  thinner  than  0.15  fim  grow 
coherently  on  GaAs  (the  mismatch  strain  is  accommodated  elastically);  between  0.15  and 
1.0  fua  the  perpendicular  lattice  parameter  (and  the  elastic  strain)  diminishes  as  misfit 
dislocations  form  to  relieve  the  strain.  Beyond  1.0  pm.  the  lattice  parameter  of  ZnSe 
matches  that  of  the  bulk  material,  so  that  the  strain  is  fully  relieved. 

Due  to  the  presence  of  the  residual  elastic  strain  in  the  films  for  h  ^  1  fim,  the 
elastic  constants  of  the  material  might  be  altered.  To  date  no  measurements  of  elastic 
constants  of  epitaxial  ZnSe  films  have  been  reported. 

In  this  paper  we  will  present  Brillouin  scattering  spectra  for  ZnSe(001)  films  grown 
by  molecular-beam  epitaxy  (MBE)  on  GaAs<001)  substrates.  The  thickness  range  of  the 
films  reflects  the  progressive  strain  relief  demonstrated  in  Ref.  1.  The  data  reveal  no 
modification  of  the  elastic  constants  with  increasing  strain;  however  the  elastic  constants 
of  the  films  are  in  good  agreement  with  those  of  bulk  ZnSe. 

Brillouin  scattering  has  been  used  successfully  to  measure  the  elastic  properties  of 
thin  films  for  a  number  of  years.2  In  thin  films,  light  is  inelastically  scattered  from 
thermally  excited  acoustic  phonons  guided  by  the  film.  Surface  and  guided  acoustic 
waves,  i.e.,  the  Rayleigh,  and  Lamb  or  so-called  Sezawa  modes,  are  localized  modes 
travelling  parallel  to  the  film  on  a  substrate.  Their  displacement  fields  show  oscillating 
amplitudes  within  the  film  and  an  exponential  decrease  of  the  amplitudes  into  the 


substrate,  essentially  confining  the  mode  energy  to  the  film.  For  large  film  thicknesses, 
the  Rayleigh  mode  is  confined  to  the  surface  of  the  film.  With  increasing  film 
thickness,  its  sound  velocity  approaches  rapidly  the  Rayleigh  sound  velocity  of  semi¬ 
infinite  ZnSe.  The  modes  can  be  described  by  partial  waves  in  each  material,  by 
solving  the  wave  equations,  and  by  taking  the  mechanical  boundary  condition  at  the 
film  surface  and  at  the  interface  between  the  film  and  the  substrate  into  account.’ 
Although  the  stress  fields  do  not  directly  enter  the  equation,  they  can  modify  the  elastic 
constants  via  higher  order  terms.  In  addition  to  the  Rayleigh  mode.  Lamb  or  so-called 
Sezawa  modes  exist  if  vyl-  <  vy$.  where  vyl*  and  vyS  are  the  bulk  transverse  sound 
velocities  of  the  layer  and  the  substrate,  respectively.  The  corresponding  sound 
velocities  vj  are  limited  to  the  regime  vyl*  <  vj  <  vy$. 

The  series  of  ZnSe  films  used  in  this  work  were  grown  on  semi-insulating 
GaAs(001)  substrates  in  a  Perkin-Elmer  Model  430  MBE  system,  as  described 
elsewhere.4-5  The  thicknesses  of  the  films  investigated  here  are  0.23.  0.53.  1.0  and  2.0 
/un,  as  measured  by  a  Dektak  surface  profiler. 

The  Brillouin  scattering  experiments  were  performed  by  measuring  surface  phonons 
propagating  both  in  [100]  and  [110]  directions  at  room  temperature  in  a  backscattering 
geometry  using  a  (3+3)-pass  tandem  Fabry-Perot  interferometer.  A  single-mode 
514.5  nm  Ar+  ion  laser  was  used  with  an  incident  power  of  up  to  300  mW,  depending 
on  the  sample  thickness  and  incident  angles.  For  sample  thicknesses  less  than  1 .0  /im  or 
an  incident  angle  less  than  50°,  the  power  was  reduced  in  order  to  avoid  surface 
damage.  The  incident  light  was  p- polarized.  The  polarization  of  the  scattered  light 
was  not  analyzed  in  order  to  achieve  a  sufficiently  high  scattering  signal.  The  aperture 
of  the  collecting  lens  (f  -  50  mm)  was  1:2.8.  Larger  aperture  did  not  allow  adequate 
resolution  of  the  localized  modes  in  the  spectra. 

Figure  1  shows  Brillouin  spectra  (full  lines)  for  a  2.0  fim  sample  (Fig.  la)  and  a 
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The  angle  of  incident  light  was  60°.  The  scattering  intensities  show  the  Stokes/anti- 
Stokes  symmetry,  which  is  characteristic  for  phonon  excitations  in  the  GHz  regime  at 
room  temperature.  The  most  intense  inelastic  peak  in  both  spectra  at  7.0  GHz  is  due  to 
the  Rayleigh  surface  mode  of  the  ZnSe  film.  The  next  intense  peak  at  13.6  GHz  in  Fig. 
la.  which  is  strongly  reduced  in  intensity  for  the  0.23  fim  sample  (Fig.  lb),  is  identified 
as  a  longitudinal  guided  acoustic  phonon  described  elsewhere.6  The  two  peaks  between 
8  and  10  GHz  in  Fig.  lb  are  the  localized  film  modes  (Lamb  modes). 

Dividing  the  measured  mode  frequencies  by  the  component  qj|  of  the  phonon  wave 
vector  parallel  to  the  surface,  which  is  determined  by  the  scattering  geometry,  one 
obtains  the  phonon  phase  or  sound  velocities.  They  are  plotted  in  Fig.  2  as  a  function 
of  q||h.  where  h  is  the  total  thickness  of  the  film.  Each  symbol  denotes  a  different 
thickness  of  the  ZnSe  film.  A  set  of  data  points  for  each  sample  was  obtained  by 
varying  q||.  i.e.,  by  changing  the  angle  of  incidence  from  4(P  to  70°  in  10°  increments. 
The  choice  of  sample  thicknesses  h  and  of  q||  allows  the  phonon  velocities  to  be 
followed  over  a  continuous  range  of  q||h.  The  error  in  the  sound  velocities  is  2%  and 
the  error  in  the  film  thickness  is  5%. 

In  order  to  determine  the  elastic  constants  of  the  film,  we  used  the  following 
procedure:  from  the  frequency  shifts  of  the  Rayleigh  mode  and  the  longitudinal  guided 
acoustic  phonon*  of  the  2.0  nm  sample  measured  in  [100]  and  [110]  directions,  we 
determined  their  corresponding  sound  velocities  and  thus  their  elastic  moduli  cr  and 
cl-  Since  the  film  thickness  is  much  larger  than  the  acoustic  wavelength  and  thus  the 
penetration  depth  of  the  Rayleigh  mode,  we  could  assume  that  the  observed  Rayleigh 
sound  velocity  is  that  of  semi-infinite  ZnSe.  The  sound  velocity  of  the  longitudinal 
guided  acoustic  phonon  is  a  very  good  approximation  of  that  of  the  corresponding 
longitudinal  bulk  phonon.*  From  this  set  of  four  experimental  values  we  determined 
the  elastic  constants  of  ZnSe  by  a  least-squares  fit.  The  residual  sum  of  squares  was 

much  smaller  than  the  experimental  error,  confirming  this  procedure.  The  results  are 
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E  listed  in  Table  1  for  comparison  with  literature  data.7*13  The  obtained  elastic  constants 

of  ZnSe  films  are  in  excellent  agreement  with  the  literature  values  obtained  for  the  bulk 

I 

materiel,  thus  ruling  out  modifications  due  to  the  Him  growth  within  5%. 

» 

The  dispersion  curves  (full  lines  in  Fig.  2)  were  calculated  using  the  elastic 
contimnim  model  described  above  with  the  obtained  elastic  constants  of  ZnSe.  Figure  2 
^  shows  the  dispersion  curves  and  the  experimental  data  in  both  [100]  and  [110] 

directions.  All  experimental  values,  especially  those  for  the  0.23  Am  sample,  are  in 
very  good  agreement  with  the  dispersion  curves  calculated  by  using  the  data  of  the  2.0 
,  Am  sample.  This  shows  that  there  is  no  change  in  the  elastic  constants  of  more  than  the 

experimental  error  between  h  -  0.23  fim  and  h  -  2.0  Am. 

In  order  to  further  test  our  results,  and  especially  to  determine  why  all  theoretically 
i  expected  modes  could  not  be  found  experimentally,  we  also  calculated  the  Brillouin 

scattering  crocs-section.  The  model  used  here  has  been  described  elsewhere13  and  takes 
both  ripple  and  elasto-optic  scattering,  as  well  as  corresponding  interference  terms,  into 
1  account.  The  photo-elastic  coefficients  have  been  taken  from  literature.11  The  sign  of 

the  coefficients  had  to  be  chosen  negative  in  order  to  produce  agreement  between  the 
calculation  and  the  experimental  data.  Figure  1  shows  the  theoretical  curves  as  dashed 
lines.  They  ere  in  good  agreement  with  the  experimental  spectra.  In  particular,  the 
relatively  small  theoretical  cross-section  for  the  first  Lamb  mode  in  the  0.S3  fim  sample 
reproduces  the  experimental  results. 

Although  transparent  films  in  general  allow  for  strong  elasto-optic  scattering 
contributions,  the  high  refractive  index  of  ZnSe  also  implies  a  large  ripple  scattering 
term.  The  dominant  scattering  contribution  for  the  Rayleigh  mode  is  provided  by  the 

> 

ripple  scattering  mechanism  at  the  film  surface.  Lamb  modes  show  also  a  significant 
contribution  from  ripple  scattering  at  the  film-substrate  interface.  For  some  modes,  e.g.. 
the  first  Lamb  mode  in  the  0.S3  fim  sample,  elasto-opticaliy  scattered  light  from  the  film 
is  comparable  in  amplitude  to  the  ripple  contributions  and  interferes  destructively  with 
light  scattered  by  the  surface  ripple,  accounting  for  the  observed  small  cross  section. 
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In  conclusion,  we  have  determined  the  elastic  constants  for  epitaxially  grown 
ZnSe(001)  films  on  GaAs(001).  No  modification  of  the  elastic  constants  could  be  found 
as  compared  to  the  literature  (bulk)  values.  In  particular,  the  elastic  constants  show  no 
measureable  thickness  dependence  for  0.23  pm  £  h  £  2.0  pm.  Thus,  the  influence  of 
the  strain  field  in  the  elastic  constants  seems  to  be  negligible.  Since  changes  in  the 
elastic  constants  are  related  to  changes  in  the  free-energy  density,  the  latter  seems  to  be 
unaffected  by  the  strain  field,  which  in  turn  is  corroborated  by  the  absence  of 
dislocation  formation  for  h  <  0.15  pm  and  the  slow  onset  for  0.15  pm  <  h  <  pm.1 
Calculations  of  the  scattering  cross  section  are  in  good  agreement  with  the  experimental 
data,  and  exhibit  scattering  intensities  mostly  dominated  by  ripple  scattering.  The  lack 
of  observation  of  some  Lamb  modes  could  be  explained  by  destructive  interferences  of 
scattering  contributions,  due  to  ripple  scattering  at  the  film  surface  and  elasto-optic 
scattering  in  the  film. 
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Table  Captions: 


Table  1.  Elastic  constants  of  ZaSe  in  comparison  to  literature  data.  All  values  are 
given  in  units  of  GPa. 

Figure  Captions: 

Figure  1.‘  Experimental  scattering  intensities  of  (a)  a  2.0  ftm  thick  sample  (b)  a  0.23  nm 
thick  sample  with  the  wave  vector  in  [100]  direction  (solid  lines)  and 
corresponding  theoretical  cross  sections  (dotted  lines).  The 
angle  of  incident  light  was  6 <f. 

Figure  2.  Dispersion  curves  in  (a)  [100]  and  (b)  [110]  directions,  respectively,  with 
experimental  data.  The  error  is  2%  for  the  sound  velocities  and  5%  for  the 
film  thicknesses  Circles:  h  -  0.23  pm.  squares:  h  -  0.53  urn.  and  triangles: 


Method  Piezo-  Ultrasonic  Ultrasonic  Brillouin  Thermal  Brillouin  Present 

resonance7  Pulse-echo*  interfere-  Scattering10  Brillouin  Resonance  Work 
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EVIDENCE  FOR  THE  EXISTENCE  OF  GUIDED  LONGITUDINAL 
ACOUSTIC  PHONONS  IN  ZnSe  FILMS  ON  GaAs 
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Abatract: 

Wa  report  theoretical  as  wall  as  experimental  evidence  for  a  new  class  of  long- 
wavelength  acoustic  film  excitations  which  wa  identify  as  longitudinal  guided  acoustic 
phonons.  Doepite  a  coupling  to  propagating  transverse  film  modes,  these  new  modes 
have  lifetimes  comparable  to  the  Rayleigh  mode  and  are  the  longitudinal  counterpart  to 
the  well-known  Lamb  or  Sezawa  modes,  with  analogous  existence  enters.  The 
properties  of  the  new  modes  are  discussed  in  detail  for  ZoSe  films  on  GaAs.  the 
material  system  for  which  we  were  able  to  observe  these  modes  by  Brillouin 
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In  the  field  of  long-wavelength  acoustic  surface  and  film  excitations,  the  Rayleigh 
mode,  which  always  exists  on  the  surface  of  a  semi-infinite  material,  is  the  most 
prominent  member.1  For  a  film  on  a  substrate,  higher-order  modes,  which  are  called 
Lamb  or  Sezawa  modes,  exist  in  the  sound  velocity  regime  bounded  by  the  transverse 
sound  velocities  of  the  film  (Vf^)  and  the  substrate  (vyS)  for  the  case  vj*  <  vj3.1  All 
of  those  modes  consist  of  shear  vertical  and  longitudinally  polarized  partial  waves 
whose  field  components  propagate  parallel  to  the  surfaces  and  decay  exponentially  with 
distance  into  the  substrate,  thus  confining  the  mode  energy  to  the  immediate  vicinity  of 
the  film  (or  surface  for  a  Rayleigh  wave).  Within  the  film,  the  shear  partial  waves  are 
propagating  in  character  and  the  longitudinal  partial  waves  are  always  evanescent. 
Recently  the  longitudinal  counterpart  to  the  Rayleigh  mode,  the  so-called  longitudinal 
resonance  with  a  souau  velocity  slightly  below  the  longitudinal  sound  velocity  and  with 
longitudinal  components  decaying  exponentially  into  the  medium,  was  predicted 
theoretically  and  obeerved  experimentally  in  gold.*"*  However,  since  this  mode  radiates 
usually  quite  strongly  into  transverss  substrate  modaa.  the  mode  is  quasi-localized,  i.e.. 
leaky,  with  a  relatively  large  liaewidth  whan  probad  by  Brillouin  spectroscopy. 
Conaaquantiy.  it  waa  expected  that  tha  longitudinal  counterparts  of  the  Lamb  (Sezawa) 
modaa.  whose  velocities  are  even  higher  than  those  of  the  longitudinal  resonances, 
would  have  even  larger  linewidthe  and  be  more  difficult  to  observe.  In  this  paper  we 
report  an  the  successful  observation  of  thaee  modes  by  Brillouin  spectroscopy  in 
epitaxial  ZaSe  films  on  GaAs  substrates.  Contrary  to  initial  expectations,  we  also 
found  that  the  lifetime  of  them  modes  can  be  long  and  the  linewidths  narrow, 
comparable  to  those  of  the  Rayleigh  wave  itself. 

The  nature  of  the  longitudinal  guided  acoustic  waves  can  be  easily  understood  by 
examining  the  acoustic  reflections  at  the  two  film  boundaries.  Since  the  principal 
contribution  to  the  total  displacement  field  consists  of  propagating,  longitudinally- 
polarized.  partial  waves  in  the  film,  which  are  incident  and  reflected  from  the  film 
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boundaries  at  grazing  angles,  mode  conversion  into  the  shear-polarized  standing  waves 
in  the  film  is  weak.  For  v\J<  v  <  v^s  where  v.  and  v^s  are  respectively  the 
film,  guided  wave,  and  substrate  longitudinal  wave  velocities,  the  longitudinal  fields 
decay  evanescently  into  the  substrate  and  no  radiation  into  substrate  longitudinal  modes 
can  occur.  Furthermore,  the  grazing  incidence  angle  condition  also  ensures  very  weak 
radiative  coupling  into  substrate  shear  modes.  These  waves  can  be  identified  as 
longitudinal  guided  acoustic  waves  because  their  properties  and  existence  conditions  are 
completely  analogous  to  those  of  the  Lamb  modes.  With  increasing  film  thickness  the 
sound  velocities  decrease  approaching  asymptotically  vyj.  Varying  the  substrate  sound 
velocity  in  a  computer  simulation  for  a  ZnSe  film  on  a  (001)  GaAs  substrate,  we  also 
found  evidence  for  higher-order  longitudinal  guided  modes. 

The  materials  system  which  satisfies  the  above  requirements  on  acoustic  properties 
conaialod  of  (001)  ZnSe  films  grown  on  semi-insulating  (001)  GaAs  substrates,  as 
described  elsewhere.*'7  Sample  thicknesses  were  0.23.  0.S3.  1.0  and  2.0  /im.  The 
Brillouin  measurements  were  performed  in  air  at  room  temperature  using  a  Sandercock- 
type  (3+3)-pass  tandem  Fabry-Perot  interferometer  in  a  beckscattering  geometry.  The 
light  source  was  a  single  mode  514.5  am  Ar+  ion  laser.  The  laser  power  at  the  sample 
was  typically  200  mW.  The  light  scattering  geometry  was  chosen  to  investigate  modes 
with  their  acoustic  wave-vector  component  parallel  to  the  surface.  q|  | .  lying  along  the 
[100]  and  [110]  directions  of  ZnSe.  The  magnitude  of  q| ,  was  varied  by  using  angles 
of  light  incidence  9  between  40*  and  70*  to  the  surface  normal.  The  aperture  of  the 
collecting  lens  was  1:2.8.  The  incident  light  was  p- polarized;  no  polarization  analysis 
of  the  scattered  light  was  made. 

Figure  1  shows  typical  spectra  for  9  -  60*  and  q|  |  in  the  [100]  direction.  The 
peak  at  7  GHz  is  identified  as  scattering  from  the  Rayleigh  phonon  (labeled  "R"  in  Fig. 
1).  The  peaks  between  about  8  and  12  GHz  are  due  to  the  localized  Lamb  (Sezawa) 
modes  (labeled  NSM).  A  pronounced  peak  of  magnitude  comparable  to  the  Rayleigh 


wave  line  is  observed  near  14  GHz  in  the  2  /mu  and  the  1  iim  films.  It  is  identified 
as  scattering  from  the  longitudinal  guided  acoustic  phonon  (labeled  MLH  in  Fig.  1).  With 
decreasing  film  thickness,  its  intensity  decreases,  accompanied  by  a  linewidth 
broadening.  -For  h  »  2pm.  the  frequency  position  of  the  mode.  14  GHz.  is  very  close 
to  the  frequency  of  a  longitudinal  acoustic  phonon  in  ZnSe  travelling  parallel  to  the 
Him.  With  decreasing  Him  thickness,  the  mode  frequency  shows  a  slight  increase. 

We  have  calculated  the  Brillouin  scattering  cross  section  using  a  theory  which  takes 
into  account  both  the  mechanical  and  the  electro-magnetic  boundary  conditions  at  the 
film  interfaces  including  both  ripple  and  elasto-optic  coupling  mechanisms.*  The  elastic 
constants  of  ZnSe  were  evaluated  from  (i),  the  frequency  shifts  of  the  peaks  associated 
with  the  longitudinal  guided  mode  which  yield  to  a  good  approximation  (see  below)  the 
sound  velocity  of  a  longitudinal  bulk  phonon  travelling  parallel  to  the  film,  and  (ii).  the 
Rayleigh  mode  in  the  2  am  film,  all  measured  in  [100]  and  [110]  directions.*  Their 
values  of  cu  ■  87.0  t  2.1.  cu  -  54.7  t  2.1  and  cM  -  39.1  ±  1.6  GPa  are  in  good 
agreement  with  values  reported  for  bulk  ZaSe.1**”  The  photo-elastic  coefficients  and 
the  dielectric  constants  of  ZnSe  and  GaAs  have  been  taken  from  literature.  14*l,*1T  In 
Fig.  1  the  calculated  cross  section  is  shown  by  dashed  lines.  Good  agreement  between 
the  theoretical  curves  and  the  experimental  data  was  achieved  by  convolving  the 
theoretical  curves  with  the  experimental  spectrometer  function.  Only  the  absolute  value 
of  the  theoretical  croas  section  has  been  adjusted.  Since  the  calculation  takes  all 
boundary  conditions  into  account,  it  is  not  surprising  that  the  new  guided  longitudinal 
phonon  is  reproduced  quite  well. 

Figure  2  shows  for  the  guided  longitudinal  mode  the  calculated  variation  of  the 
squared  parallel  and  perpendicular  components  of  the  displacement  field  in  the  film 
(negative  abscissa)  and  the  substrate  (positive  abscissa)  normalized  with  respect  to  the 
maviwMim  value  of  the  parallel  component  in  the  substrate.  They  are  a  measure  of  the 
displacement  energy.  In  the  film  the  displacements  are  dominantly  longitudinal  and 


their  squared  values  exceed  those  of  the  substrate  by  a  factor  of  70  for  h  -  2  (tm.  The 


two  pairs  of  partial  waves  in  the  film,  longitudinal  and  shear,  each  form  standing  wave 


patterns.  The  larger  periodicity  is  due  to  the  longitudinal  partial  waves,  whereas  the 


shorter  n«ein*tirtM  are  due  to  the  transverse  components.  The  inset  of  Fig.  2  indicates 


the  propegation.chrsctions  associated  with  the  longitudinal  (L)  and  transverse  CD  Him 


and  substrate  partial  waves.  The  longitudinal  partial  waves  travel  nearly  parallel  to  the 


film  whereas  the  transverse  ones  have  a  larger  wave  vector  component  perpendicular  to 


the  film. 


As  an  existence  criterion  far  guided  waves,  it  is  necessary  that  the  partial  waves 


interfere  constructively.  Thus  the  wave-vector  component  perpendicular  to  the  Him. 


qj_.  must  obey  tho-resonance  condition 


2  qj_  h  *  *  5I  ■  2  n  r. 


where  5$.  Jj-  are  the  phase  shifts  upon  reflection  at  the  two  film  surfaces  respectively. 


and- the  mode  index  n  is  a  positive  integer.  A  calculation  of  the  phase  shifts  using  a 


simplified  model1*  (Le.  elastically  isotropic  medium)  shows  that  Eq.  1  is  fulfilled  very 


accurately  for  the  longitudinal  partial  waves  with  n  -  1.  Eq.  1  also  predicts  several 


additional  properties  of  the  mode.  Because  the  elastic  properties  of  the  media  bounding 


the  film  are  dissimilar,  there  is  a  minimum  film  thickne 


below  which  Eq.  1 


cannot  be  satisfied.  With  increasing  film  thickness  the  value  of  qj_  decreases  in  order  to 
fulfill  Eq.  1.  thus  decreasing  the  modal  sound  velocity.  A  decreasing  value  of 
implies  a  decreasing  amount  of  mode  conversion  into  shear  waves  resulting  in 


increasing  localization,  as  experimentally  observed  (sec  Fig.  1  ).  In  the  limit  of  infinitely 


thick  films.  qj_  approaches  zero,  which  implies  that  the  guided  longitudinal  acoustic 


phonon  becomes  an  ordinary  longitudinal  bulk  phonon  travelling  parallel  to  the  film.  It 


should  be  stressed  that,  for  a  film  of  finite  thickness,  no  longitudinal  bulk  phonon 
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travelling  parallel  to  the  film  may  exist  because  the  mechanical  boundary  conditions  at  ! 

i 

the  interfaces  cannot  be  satisfied. 

li 

Fig.  3  shows  results  of  computer  simulations  of  the  maximum  of  the  squared  * 

parallel  displacement-field  amplitude  of  the  longitudinal  guided  mode,  normalized  as  in  I 

Fig.  Z  as  a  function  of  the  film  thickness.  With  increasing  film  thickness  the  sound 
velocity  shows  an  over-all  decrease  with  small  superimposed  oscillations.  The  intensity 

of  the  mode  increases,  accompanied  by  strong  aperiodic  oscillations  which  exhibit  j 

i 

monotonicaUy  increasing  maxima  with  intervening  regions  of  smaller  maxima,  j 

accompanied  by  corresponding  oscillations  in  the  linewidth.  The  oscillatory  behavior  is  ! 

due  to  the  mode  conversion  at  each  reflection  of  transverse  into  longitudinal  partial  ! 

waves  and  vice  versa.  The  transverse  waves  can  also  fulfill  a  resonance  condition 

I 

similar  to  Eq.  1.  however  with  a  large  n  which  also  increases  with  increasing  film 

i 

thickness.  Since  the  total  transmission  of  energy  across  the  interface  is  proportional  to 
the  transmission  coefficient  at  the  interface  and  the  displacement  amplitude  of  the 
transverse  film  partial  waves,  a  film  resonance  of  the  transverse  modes  implies  a 
maximum  in  the  energy  loss  into  the  substrate.  Since  the  transverse  mode  resonance 
condition  is  periodically  fulfilled  with  increasing  film  thickness,  the  guided  mode 
attenuation  also  oscillates,  accompanied  by  a  corresponding  oscillation  in  the  linewidth. 

The  two  different  types  of  maxima  can  be  related  to  the  odd  or  even  value  of  n:  a 
calculation  in  the  iaotropic  approximation1*  of  the  phase  factors  for  the  reflection  of  the 
transverse  film  inodes  at  the  surface  (&j  -  -v/2)  and  at  the  interface  (ft;  a  0.2)  yields  a 
major  maximum  in  the  displacement  Held  for  even  n  and  a  secondary  maximum  for  odd 
a. 

In  summary,  we  have  found  a  longitudinal  guided  acoustic  mode  in  ZoSe(00l)  films 
on  GaAs(00l).  Its  calculated  linewidth  is  unusually  small  typically  0.5%  of  its  mode 
frequency.  Despite  the  decay  channel  into  transverse  substrate  modes,  the  new  mode  is 
the  longitudinal  counterpart  to  the  (transverse)  Lamb  modes  in  a  film.  Thus  the 
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classification  of  long- wavelength  acoustic  surface  and  film  excitations  now  exhibits 
complete  symmetry  with  respect  to  shear  and  longitudinal  waves,  and  includes  surface 
modes  (Rayleigh  mode,  longitudinal  resonance)  and  film  modes  (Lamb  modes,  guided 
longitudinal  modes).  The  new  longitudinal  guided  wave  reported  here  offers  a  number 
of  interesting  possibilities.  For  example,  it  allows  a  very  accurate  determination  of  the 
longitudinal  sound  velocity  in  films,  a  parameter  which  is  otherwise  difficult  to 


evaluate.  Here  we  actually  utilized  this  feature  to  evaluate  the  elastic  constants  of  the 
ZnSe  film.  Since  the  displacement  field  is  nearly  completely  polarized  parallel  to  the 
film,  the  mode  can  be  used  to  monitor  effects  which  "wviniat*  only  dilatational  waves, 
and  vice  versa.  The  observation  and  application  of  the  novel  mode  in  systems  other 
than  ZnSe/CaAs  is  presumably  only  restricted  by  the  existence  criteria  on  the  acoustic 
velocities  described  above. 
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FIOURE  CAPTIONS 


Brillnuin  spectra  of  ZnSa(001)  films  on  GaAs(001)  substrate.  The  angle  of  incidence  is 
6(P.  the  fihn  thicknesses  are  as  indicated.  R:  Rayleigh  mode.  S:  Lamb  (Sezawa)  modes. 
L  icagiinjinal  ytyiarf  niode 


Ql * 

‘Modal  dispiecement  Held  profile  for  a  2  fan  ZaSe(00l)  film  on  GaAs(001),  normalized 
with  respect  to  the  maxfmmn  value  of  the  parallel  component  in  the  substrate.  The 
aepdve  values  of  the  ahsciiM  refer  to  the  film,  whereas  the  positive  values  refer  to  the 
substrate.  The  ordinate  is  the  squared  amplitude  component  parallel  to  the  Him  (full 
lines)  and  perpendicular  to  the  film  (broken  lines)  of  the  displacement  field.  Inset: 
Sketch  of  the  pertial  waves  which  constitute  the  longitudinal  guided  mode.  T: 
transverse.  L:  longitudinal  partial  waves.  Not  shown  is  the  exponentially  decaying 
substrata  longitudinal  partial  wave. 


SfcJi 

Surface  plot  and  contour  plot  of  the  squared  parallel  displacement  field  as  a  function  of 
the  film  thickness  h  and  the  sound  velocity  v. 
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BRILLOUIN  SCATTERING  FROM  COLLECTIVE  SPIN  WAVES 
IN  MAGNETIC  SUPERLATTICES  (invited) 
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A.  Boufelfel,  C.M.  Falco 

Dapt.  of  Physics  and  Optical  Sciences  Center.  University  of  Arizona.  Tucson.  AZ  85721 
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E.  Zimgiebl,  J.D.  Thompson 

Los  Alamos  National  Laboratory.  Los  Alamos,  NM  87545 


We  report  on  the  observation  and  the  analysis  of  collective  magnetostatic  spin-wave 
excitations  in  magnetic  superlattices.  The  influence  of  interface  anisotropies,  which  can 
become  dominant  for  small  modulation  wavelengths,  is  discussed.  For  the  system  Fe/Pd 
we  show  that  Brillouin  spectroscopy  experiments  in  combination  with  the  measurement 
of  the  saturation  magnetization  by  a  SQUID  magnetometer  give  evidence  for  a  magnetic 
polarization  of  the  Pd  spacer  layers,  as  well  as  for  a  small  negative  out-of-plane 
interface  anisotropy  constant  of  K,  -  -0.15  erg/cm3. 


PACS-Numbers:  75.30.Ds,  75.30.Gw,  75.70.-i.  87.35,+c 


Magnetic  multilayered  structures  have  gained  considerable  interest  in  the  past  few  years 
due  to  their  unusual  magnetic  properties.  Besides  novel  static  phenomena  like  the  anti¬ 
phase  domains  in  Gd/Y  superlattices1  and  the  spin  spiral  effects  in  Dy/Y  superlattices2 
the  so-called  collective  spin  wave  excitations  are  the  first  and  the  only  dynamic 
collective  features  of  superlattices  reported  so  far. 2-10  In  most  cases  they  are  made  up 
from  surface  spin-waves  n  each  magnetic  layer  coupled  by  their  dipolar  stray  Helds 
via  the  intervening  nonmagnetic  spacer  layers  giving  rise  to  a  new  kind  of  spin-wave 
band,  which  exhibits  partially  bulk-mode-  as  well  as  surface-mode-like  behavior.  A 
study  of  these  excitations  can  reveal  information  about  the  magnitude  and  the  direction 
of  the  saturation  magnetization,  about  magnetic  anisotropies,  and  about  a  possible 
polarization  of  the  spacer  material.  In  particular  Brillouin  spectroscopy  seems  to  be  best 
suited  to  study  magnetic  interface  anisotropies  due  to  the  sensitivity  of  spin  waves  to 
the  anisotropy  Helds.  It  is  therefore  of  crucial  importance  to  consider  Brillouin 
scattering  experiments  in  addition  to  established  methods  of  magnetic  measurements  in 
order  to  gain  more  information  about  the  magnetism  of  multilayered  structures. 

So  far  there  are  only  few  reports  on  the  observation  of  spin  waves  in  magnetic 
super  lattices  by  Brillouin  spectroscopy.6-10  After  theoretical  predictions  of  the  existence 
of  the  collective  modes  by  Camley  et  al..2  Gnlnberg  and  Mika.4  and  by  Emtage  and 
Daniel.2  first  experimental  evidence  was  provided  by  Brillouin  scattering  experiments  in 
Mo/Ni  superlattices  by  Grimsditch  et  al.6  and  Kueny  et  al.7.  A  first  detailed 
experimental  proof  of  the  predicted  magnetic  properties  of  the  collective  spin-wave 
band  including  measurements  of  the  band  shapes  was  presented  by  Hillebrands  et  al. 
for  Fe/Pd  and  Fe/W  superlattices.  *•’ 

In  this  paper  we  present  Brillouin  scattering  experiments  on  Fe/Pd  superlattices  in 
addition  to  a  brief  review  of  earlier  reported  results.*’9  We  outline  how  to  fully  include 


interface  anisotropies  in  the  calculation  of  the  spin-wave  frequencies  and  re-analyze 
within  the  new  model  the  experimental  data'  of  the  Fe/Pd-system.  We  demonstrate  that 
a  previous  description**9  of  the  observed  spin- wave  frequencies  by  a  decrease  of  the 
saturation  magnetization  4SM,  for  small  modulation  wavelengths  A  can  be  improved  by 
the  approach  of  considering  a  small  out-of-plane  interface  anisotropy  assuming  a  value 
of  4ffM,  -  17  ±  2  kO  independent  of  A.  A  recently  reported  strong  increase  of  4*^ 
for  small  A.  as  measured  by  den  Broader  et  al.  using  a  vibrating  sample  magnetometer,11 
can  consistently  be  described  by  the  assumption  of  a  magnetic  polarization  of  the  Pd 
spacer  layers.  The  interpretation  is  corroborated  by  magnetization  measurements  in 
magnetic  fields  parallel  and  perpendicular  to  the  layers  using  a  SQUID  magnetometer. 

The  Brillouin  scattering  experiments  were  performed  at  room  temperature  under 
vacuum  in  backscattering  geometry  using  a  (3+3)- pass  tandem  Fabry-Perot 
interferometer  described  elsewhere.1*  A  single-mode  S14S  A  Ar+  ion  laser  was  used 
with  an  incident  power  of  p-polarized  light  of  up  to  250  mW.  The  sampling  time  per 
spectrum  was  typically  5  h.  The  inelastically  scattered  light  was  analyzed  depolarized  in 
order  to  suppress  surface  phonon  signals.  The  applied  magnetic  field  was  oriented 
parallel  to  the  layer  planes  and  perpendicular  to  the  scattering  plane. 

The  Fe/Pd  superlattice  samples  were  prepared  on  single  crystal  sapphire  substrates 

using  a  sputtering  technique  described  elsewhere.14  The  sample  parameters  are  listed  in 

0 

Ref.  8.  The  modulation  wavelength  A  varies  between  33  and  188  A,  and  the  number 

of  bilayers  N  varies  between  49  and  90.  The  modulation  wavelengths  of  the  different 

samples  were  determined  by  x-ray  diffraction.  As  shown  by  Bragg  and  wide-film 

Debye-Scherrer  x-ray  diffraction  the  layers  grew  with  a  preferred  orientation  of  bcc 

Fe(l  10)  planes  and  fee  Pd(lll)  planes,  with  no  in-plane  orientation.  The  samples 

0 

exhibited  long-range  structural  coherence  of  at  least  300  A  perpendicular  to  the  layers. 
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In  Figs.  1  and  2  typical  Brillouin  spectra  of  magnetic  superlattices  are  displayed  for  an 

applied  magnetic  field  of  1  kC.  Fig.  la  shows  the  spectrum  of  a  Fe/Pd  superlattice  with 

0 

a  modulation  wavelength  of  A  -  46.2  A.  The  thickness  d  of  the  magnetic  material  is 
close  to  that  of  the  nonmagnetic  material  (dj.  The  band  of  collective  spin-wave 
excitations  can  clearly  be  identified  in  the  right-hand  part  of  the  spectrum  by  its 
specific  asymmetric  shape.  The  density  of  states  is  largest  at  small  frequency  shifts  and 
decreases  asymmetrically  toward  the  upper  band  edge.  At  the  latter  discrete  spin-wave 
modes  can  still  be  resolved  due  to  the  small  thickness  and  the  finite  number  of 
bilayers.***  Their  strong  Stokes/anti-Stokes  asymmetry  identifies  them  as  surface-mode- 
like  spin  waves.  On  the  other  hand,  the  modes  near  the  lower  edge  of  the  spin-wave 
band  are  bulk-mode-like  because  of  their  small  Stokes/anti-Stokes  asymmetry.  If  we 
neglect  the  discrete  modes  near  the  upper  band  edge  the  shape  of  the  spin-wave 
excitation  band  is  qualitatively  similar  to  the  calculated  Brillouin  scattering  cross  section 
for  the  semi-infinite  superlattice  system  Mo/Ni.* 

0 

In  Fig.  lb  we  show  the  Brillouin  spectrum  of  an  Fe/Pd  superlattice  with  d  -  41.7  A 

and  a  much  larger  do  -  138.7  A.  In  this  case  the  spin-wave  band  becomes  narrower 

and  no  discrete  peaks  are  observable  near  the  upper  band  edge.  A  drastically  different 

0  0 

spectrum  is  found  for  the  case  of  do  (9.1  A)  much  smaller  than  d  (41.0  A)  as  shown  in 
Fig.  Ic.  In  this  case  a  very  intense  discrete  mode  is  found  near  27.7  GHz  in  the  anti- 
Stokes  spectrum.  This  superlattice  surface  spin-wave  mode,  which  travels  about  the 
total  superiattice  stack,  is  allowed  to  exist  besides  the  collective  spin-wave  band. 1-5  It 
would  merge  with  the  latter  for  d  d,. 


The  effect  of  inverting  the  applied  magnetic  field  on  the  Brillouin  cross  section  is 
demonstrated  in  Fig.  Id.  Since  the  direction  of  the  applied  field  defines  the  sense  of 
revolution  of  each  surface  spin-wave  mode  about  each  magnetic  layer,  an  inverted  field 

102 


exchanges  the  Stokes  and  anti-Stokes  parts  of  the  spectra. 

0 

Fig.  2  shows  Brillouin  spectra  of  a  Fe/Pd  superlattice  sample  with  d  -  89.4  A  and  do  - 
99.0  A  for  different  applied  magnetic  fields  H.  With  increasing  H  the  spin- wave 
frequencies  increase,  accompanied  by  a  band-width  narrowing. 

Since  the  thickness  of  each  magnetic  layer  in  the  superlattices  considered  here  is  as  low 
0 

as  16  A  for  the  smallest  modulation  wavelength,  an  influence  of  interface  anisotropies 

on  the  spin-wave  frequencies  cannot  generally  be  excluded.  For  a  single  epitaxial  layer 

of  epitaxial  Fe(110)  on  a  W(110)  substrate  prepared  in  UHV  by  electron  beam 

evaporation  and  in  situ  measured  by  Brillouin  scattering  an  increase  of  the  spin-wave 

0 

frequencies  toward  smaller  Fe  thicknesses  has  been  observed  for  d  <  60  A.15  Using  a 
modified  Damon-Eshbach  theory,  which  includes  properly  the  magnetic  in-plane  and 
out-of-plane  anisotropy  contributions,  a  quantitative  description  of  the  observed  mode 
behavior  could  be  achieved.  From  a  fit  of  the  experimental  data  the  anisotropy 
constants  could  be  determined.19 

In  the  following  we  will  describe  the  theory  used  in  the  present  work  on  magnetic 
superlattices,  which  includes  magnetic  interface  anisotropy  contributions.  It  is  an 
extension  of  a  theoretical  approach  described  elsewhere,13  which  allows  in  the 
magnetostatic  limit  for  the  calculation  of  spin- wave  frequencies  of  superlattices  under 
the  inclusion  of  volume  anisotropy  contributions  described  by  the  volume  anisotropy 
constant  K,  and  an  uniaxial  anisotropy  constant  Ku.  We  dropped  contributions  from  the 
latter,  since  there  is  no  evidence  in  the  superlattices  for  uniaxial  contributions,  which 
cannot  be  described  more  adequately  by  surface  anisotropies.  Exchange  contributions 
are  neglected.3"5*10*13 
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Since  the  Fe  layers  investigated  here  all  have  the  (1 10)  planes  parallel  to  the  interfaces 


of  the  superlattice,  we  restrict  ourselves  in  the  following  to  this  case.  We  will  use  a 


spherical  coordinate  system,  where  9  is  the  azimuthal  angle  measured  against  the 


interface  normal  and  ♦  is  the  polar  angle  measured  with  respect  to  the  crystallographic 


[001]  axis  in  the  layer  plane.  The  free  energy  density1*-1*  of  the  system  consists  of  i) 


the  potential  energy  density  fp  of  the  magnetization  Ms  in  the  external  field  H. 


f  p  -  -  H  sin  9  cos  -  <&) 


where  is  the  angle  between  the  crystallographic  [001]  direction  in  the  surface  and 


the  external  field  H  ii)  the  free  energy  density  fv  due  to  volume  anisotropies. 


^ 1  co «*  9  +  sin*  9  [  siri4  ♦  +  sin*  (2*)]  +  sin*  (29)  jcos*  $  -  ~  sin*  *  j  - ,  (2) 


and  iii)  the  free  energy  density  fj  due  to  interface  anisotropies  in  the  homogeneous 


magnetization  approximation,1’ 


fj  -  ~  (Kj  cos*  9  +  K^p  sin*  9  cos*  $)  . 


where  d  is  the  thickness  of  the  magnetic  layers  and  the  factor  2  counts  the  two 


interfaces.  fCj  and  K^  p  are  the  out-of-piane  and  the  in-plane  anisotropy  constants. 


respectively.  Since  the  films  considered  here  are  isotropic  in  plane,  we  have  to  average 


Eqs.  1-3  over  the  in-plane  angle 


r\\v% 


ij 


fp»°  -  -  Ms  H  sin  0  . 


(4) 


0  +  COS3  0  . 


iso  _  f  l  -:_4 

'  •  T  |l  sm 


fi®  -  I  [k,  erf  e  +  i  ^  sta>  aj . 


With  the  resulting  free  energy  density  f. 


f  "  fp  ^  fy  ^  f[  t 


we  are  now  able  to  determine  the  components  of  the  dynamic  susceptibility  tensor 


4  w  7*  f$$ 
Q1  sin3  0  ’ 


4  w  f  7*  . 


1  n  •  :  ..  w 

w  iirr',1'Ms- 


4  »  t3  f  00 
X22  -  q3 
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Q1  -  Wo2  -  w*  ,  (12) 


“•  •  Kgis'g  Jfee  f**  -  fe»* .  <13> 

where  f ^  refer  to  the  second  derivatives  of  f  with  respect  to  /z  and  i>.  y  is  the 
gyromagnetic  ratio  and  u>  the  spin-wave  frequency.  Since  for  all  superlattices 
considered  in  this  work  the  magnetization  has  no  out-of-plane  component  (see  below) 
we  may  assume  in  the  following  9  -  90°.  The  evaluated  expressions  for  f@@,  f<jxj>  and 
f@4  are 


fee-MjH-  §k.*  4  [*,- 

(14) 

f$$  -  Mj  H  , 

(15) 

fe*  -  0  • 

(16) 

It  should  be  pointed  out,  that  due  to  the  (110)  orientation  of  the  crystallites  in  each 
magnetic  film  the  anisotropy  contributions  do  not  average  to  zero.  However,  only  a 
linear  combination  of  the  three  anisotropy  constants  enters  the  calculation,  contrary  to 
results  of  epitaxial  films,  where  the  constants  enter  independently  in  f@@  and  f$$.ls 


The  rest  of  the  calculation  is  straightforward.  Solving  Maxwell's  equations  in  the 
magnetic  and  nonmagnetic  layers  and  taking  the  magnetic  boundary  conditions  at  each 
interface  into  account,  the  problem  of  finding  solutions  for  propagating  spin  waves  can 
be  reduced  to  finding  the  zeros  of  the  boundary-condition  determinant.  Since  at  each 


Vi? 
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interface  only  magnetic  interactions  between  the  adjacent  layers  have  to  be  considered 
the  boundary-condition  determinant  can  be  written-  hr  the- form  of  a  band  matrix.  This 
allows  to  apply  very  effective  numerical  tools  in  order  to  efficiently  calculate  the  spin- 
wave  frequencies  even  for  90  bilayers. 

The  parallel  component  of  the  spin-wave  wave-vector,  q||.  is  defined  by  the  scattering 
geometry  of  the  Brillouin  scattering  experiment.  The  perpendicular  component  qj_  is 
given  in  each  layer  by  the  relation13 


qL"  qll 


1  +  Xn  «ina  « 
1  +  Xn 


(17) 


where  a  is  the  angle  between  q||  and  4*Mg.  In  the  absense  of  surface  anisotropies  q]_is 
always  imaginary,  resulting  in  surface  spin-wave  modes.  However,  in  particular  for 
small  magnetic  thicknesses  and  large  positive  surface  anisotropy  contributions.  qj_can 
become  real  and  the  spin  waves  have  bulk-like  character  in  the  corresponding  layer. 

Figs.  3  and  4  show  results  of  the  calculations.  We  used  the  parameters  obtained  from 
the  fit  of  the  spin-wave  frequencies  in  single  epitaxial  Fe(l  10)  films  on  W(1 10)  reported 
in  Ref.  15.  The  parameters  are  given  in  Figs.  3  and  4.  Fig.  3  shows  the  obtained  spin- 
wave  frequencies  for  a  N-l 0  superlattice  as  a  function  of  the  Fe  layer  thickness  d. 
which  is  equal  to  the  thickness  do  of  the  nonmagnetic  spacer  material.  As  in  the  case  of 
single  layers  the  spin-wave  frequencies  increase  for  small  d.  The  mode  of  highest 
frequency,  which  is  the  surface  spin-wave  mode  revolving  about  the  total  superlattice 
stack,  shows  an  overall  frequency  increase  with  decreasing  d  and  is  well  separated 
from  the  residual  modes.  Contrary  to  the  case  of  vanishing  interface  anisotropies,  this 
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mode  cannot  be  described  as  a  surface  spin  wave  of  a  magnetic  layer  of  a  thickness 
equal  to  the  total  thickness  of-  the  superlattice  stack,1-5  -since  in  -the  -latter  -case  the 
surface  anisotropy  contribution  should  be  negligible  due  to  the  1/d  factor.  All  residual 
modes  are  nearly  degenerate  in  frequency  for  very  small  d.  An  analysis  of  the  wave 
solutions  in  each  layer  shows,  that  due  to  the  dominating  interface  anisotropy  terms  the 
spin  waves  are  bulk-mode-like  in  each  layer,  which  compared  to  surface  modes  have 
smaller  dipolar  stray  fields  in  the  nonmagnetic  spacer  layers  and  are  thus  less  coupled 
to  each  other. 

0 

Pig.  4  shows  the  spin-wave  frequencies  of  a  NU10  superlattice  with  d  -  d®  -  30A  as  a 
function  of  the  surface  anisotropy  constant  Kg.  In  the  plotted  range  the  energy 
contribution  of  the  surface  anisotropy  is  too  small  to  turn  the  magnetization  out  of  plane 
for  negative  values  of  Kg.  With  increasing  Kg  the  spin-wave  frequencies  increase, 
whereas  for  negative  values  of  Kg  the  spin-wave  frequencies  are  reduced. 

In  order  to  test  the  influence  of  interface  anisotropies  in  Fe/Pd  superlattices  as  well  as 
to  address  the  above  described  discrepancies  in  the  reported  values  of  the  saturation 
magnetization.''*'11  we  have  re-analyzed  the  previously  measured  Brillouin  scattering 
data  using  the  new  model.  In  order  to  support  the  reported  increased  values  of  the 
saturation  magnetization  for  vapor-deposited  Fe/Pd  superlattices  of  Fe  thicknesses 
between  2  and  6  A  (Ref.  11)  we  also  measured  for  some  of  our  samples  4?rMs  and  the 
saturation  field  (kink  field)  Hg  by  means  of  a  SQUID  magnetometer  for  applied  magnetic 
fields  H  up  to  54  kG  parallel  as  well  as  perpendicular  to  the  layers.  The  results  from 
the  SQUID  measurements  are  listed  in  Table  1  together  with  those  from  Ref.  11.  For  all 
samples  we  obtained  higher  values  of  the  saturation  magnetization  as  compared  to  the 
Brillouin  scattering  results  reported  in  Refs.  8  and  9.  Fig.  5  shows  a  plot  of  the 
saturation  magnetization  obtained  from  the  SQUID  measurements  as  a  function  of  the 


inverse  Fe  thickness  i/d  including  the  data  from  Ref.  11,  except  sample  1/8.2.  The 
values  of  4»Mj  of  each  superlattice  sample  -have  been  obtained  by  dividing  the 
measured  magnetic  moment  by  the  total  Fe  volume.  The  data  points  lie  roughly  on  a 
straight  line,  which  can  be  described  by 


4rMs<d)  .  4»Ms(d-oo) 


(18) 


where  d,  is  the  amount,  by  which  the  thickness  of  the  magnetic  layers  differs  from  that 
of  the  Fe  layers.  4»Ms(oo)  is  the  intrinsic  magnetization  of  Fe.  An  extrapolation  of  this 
line  into  the  regime  of  the  data  of  Ref.  11  agrees  with  the  latter,  in  particular  with  the 
sample  of  largest  Fe  thickness  (d  -  6  A).  The  line  is  given  by  4rMs(oo)  -  14.5  ±  1.5  kG 
and  dj  •  -1.5  A.  assuming  a  uniform  magnetization  across  the  total  magnetic  thickness. 
The  negative  value  of  d,  implies  that  the  thickness  of  the  magnetically  active  material 
is  larger  than  the  Fe  thickness,  i.e.  that  Pd  is  magnetically  polarized.  However,  since  in 
the  Pd  material  only  the  product  of  4rMs(oo)  of  Pd  and  d$  enters  Eq.  (18).  we  cannot 
determine  each  of  both  factors  independently. 

An  intrinsic  saturation  magnetization,  which  is  independent  of  A,  allows  for  a  consistent 
description  of  the  Brillouin  scattering  data  as  well.  An  increased  value  of  4irMg  would 
result  in  increased  spin- wave  frequencies,  which,  however,  are  not  observed 
experimentally.  On  the  other  hand,  a  small  increase  in  d  with  a  corresponding  decrease 
in  d«  modifies  the  spin-wave  frequencies  only  slightly. 


Assuming  therefore  the  saturation  magnetization  to  be  independent  of  the  modulation 
wavelength  A,  magnetic  anisotropy  effects  have  to  account  for  the  observed  decrease  of 
the  spin-wave  frequencies  with  decreasing  A.  Since  only  a  linear  combination  of  all 
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aaisotropy  constants  enters  the  formula  (see  Eq.  14  and  IS)  and  K^p  is  usually  much 
smaller  than  K^,  we  assumed  K^p  -  0  and  for  Kt  a  literature  value  of  4.8- 10*  erg/cm3.20 


A  good  fit  of  the  experimental  data  could  be  achieved  within  the  experimental  error  by 
a  value  of  4*Mg  -  17  *  2  kO  and  an  out-of-plane  anisotropy  constant  -  -0.15  ±  0.03 
erg/cm3.  For  comparison,  an  analysis  of  the  saturation  field  Hg  (see  Table  1)  as  a 
function  of  1/d  as  determined  by  the  SQUID  magnetometer  yields  a  value  of  -  -0.38 
±0.1  erg/cm3.  The  latter  value  has  the  same  sign  as  the  value  obtained  from  the 
Brillouin  spectroscopy  data.  The  higher  absolute  value  might  be  due  to  the  contribution 
of  higher  order  anisotropy  constants,  which  may  show  no  influence  on  the  spin-wave 
frequencies.31  However,  both  values  are  too  small  to  rotate  the  direction  of 
magnetization  out  of  plane.  This  is  confirmed  by  the  absence  of  hysteresis  effects  in  the 
perpendicular-field  SQUID  measurements. 

In  conclusion,  we  have  shown  that  the  observation  of  spin-wave  spectra  in  superlattices 
by  Brillouin  spectroscopy  reveals  important  information  about  the  magnetic  properties  of 
these  artificially  layered  materials.  In  particular  we  have  demonstrated  the  influence  of 
interface  anisotropies  on  the  spin-wave  frequencies  which  may  become  dominant  for 
small  modulation  wavelengths.  Since  the  spin-wave  frequencies  are  sensitive  to  the 
intrinsic  saturation  magnetization,  the  combination  of  Brillouin  spectroscopy  with 
magnetometric  investigations  allows  potentially  for  the  determination  of  dead  layers  as 
well  as  polarization  of  the  non-magnetic  spacer  layers  at  the  interfaces.  We  applied 
both  types  of  methods  to  the  Fe/Pd  superlattice  system.  From  the  comparison  of  the 
results  a  consistent  description  could  be  found  by  assuming  a  saturation  magnetization 
of  17  ±  2  kO,  which  is  independent  of  the  modulation  wavelength,  a  polarization  of  the 
Pd  layer  at  each  interface  and  a  small  negative  out-of-plane  interface  anisotropy  of  Kj 
-  -0.15  erg/cm3.  The  obtained  value  for  the  saturation  magnetization  is  about  20% 
smaller  than  the  bulk  Fe  value  of  21  kG,  which  we  attribute  to  the  sputtering  process. 
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Future  investigation  on  both  sputtered  and  MBE-grown  Fe/Pd  systems  will  help  to 
clarify  this  point.  Presently  the  physical  origin  of  the  Pd  polarization  is  not  clear.  The 
situation  is  reminiscent  of  the  formation  of  giant  moments,21  but  the  physical  conditions, 
i.e.  the  Fe  concentration  and  the  temperature  range,  are  not  comparable.  From  the 
Bethe-Slater  diagram  and  the  lattice  constant  mismatch,  which  causes  a  stretching  of  the 
spacing  of  nearest  neighbor  Pd  atoms  by  10  %,  the  interface  mismatch  would  not  give 
rise  to  a  magnetic  moment  in  the  Pd  layer.  On  the  other  hand,  a  hybridization  of  the 
Pd(d)  states  with  the  Fe(d)  states  is  very  likely  due  to  the  large  density  of  Pd  states  at 
the  Fermi  level,  and  might  account  for  the  observed  behavior.  A  fully  self-consistent 
band  structure  calculation  would  be  most  desirable  for  this  case. 
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TABLE  CAPTION 


Table  I:_ 

Parameters  and  results  of  the  Fe/Pd  superlattices  investigated  by  SQUID  measurements. 
npe  (np^):  number  of  atomic  Fe  (Pd)  layers  per  bilayer;  A:  modulation  wavelength; 
4aMg:  measured  SQUID  values  of  the  saturation  magnetization;  perpendicular 

saturation  Held. 


FIGURE  CAPTIONS 


Figure  1: 

Room  temperature  Brillouin  spectra  of  Fe/Pd  superlattices  in  an  applied  magnetic  field 
of  1  kG:  (a)  d  -  21.9  k  and  d,  -  24.3  k.  (b)  d  -  41.7  k  and  d,  -  138.7  A.  and  (c.d)  d 
-  41.0  k  and  d«  -  9.1  k.  In  (d)  the  direction  of  the  applied  field  has  been  reversed 
compared  to  (c).  The  angle  of  incidence  of  the  laser  beam  is  45°. 

Figure  2: 

Room  temperature  Brillouin  spectra  of  a  Fe/Pd  superlattice  consisting  of  49  bilayers 
0  0 

with  d  -  89.4  A  and  d,  -  99.0  A  for  different  applied  magnetic  fields  as  indicated  in 
the  Figure.  The  angle  of  incidence  of  the  laser  beam  is  45°. 

Figure  3: 

Calculated  spin-wave  frequencies  as  a  function  of  the  layer  thickness  for  a  N-10 
superlattice  with  equal  thicknesses  of  the  magnetic  and  nonmagnetic  layers.  The 
parameters  are  as  indicated  in  the  Figure. 
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Figure  4: 

Calculated  spin-wave  frequencies  as  a  function  of  the  out-of -plane  interface  anisotropy 
constant  K,.  The  parameters  are  as  indicated  in  the  Figure. 

Figure  5: 

Values  of  the  saturation  magnetization  obtained  from  the  SQUID  measurements  (crosses) 
and  values  from  Ref.  11  (circles)  as  a  function  of  the  inverse  Fe  thickness.  The  solid 
line  is  a  fit  to  our  data  and  the  dashed  line  is  an  extrapolation  into  the  regime  of  the 
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Table  1: 


npe/npd 

number  of 

bilayers 

A 

(A> 

4»M$ 

(kG) 

^oHk 

(kG) 

44/44 

49 

188.4 

15.6 

20.0 

20/20 

70 

86.8 

16.9 

17.2 

10/10 

90 

42.2 

19.2 

14.7 

3/8.21* 

125 

26.5 

27.0 

11.0 

2/8.2t 

136 

24.4 

27.4 

8.5 

1/8.2+ 

150 

22.2 

25.2 

4.6 

"*"  data  from  den  Breeder  et  al.11 
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FIGURE  A 


HIGH  ENERGY  3.8  MEV  BACK- 
SCATTERING  SPECTRUM  OF  A 

THIN  YBaCuO  FILM  (SAMPLE  B) 

ON  AN  ALUMINUM  OXIDE  SUBSTRATE. 
THE  Cu.Y  AND  Ba  SIGNALS  ARE 
WELL  RESOLVED.  THE  OXYGEN 
SIGNAL  IS  UNOBSERVABLE. 


TABLE  2 


RESULTS  FOR  SAMPLE  B 
ELEMENT  Nt  REL.  ATOMIC  % 

Y  2.10+0.06  23.0+0.1 

Ba  2.63+0.08  28.8+0.1 

Cu  4.40+0.13  48.2+0.2 

Nt  =  AREAL  DENSITY  (El 7  ATOMS/SQ  CM) 


FIGURE  5 


PORTION  OF  THE  HIGH  ENERGY 

4.7  MEV  BACKSCATTERING  SPECTRUM 


OF  A  THIN  YBaCuO  FILM  (SAMPLE  D) 
SHOWING  INTERFERENCE  OF  THE  Cu, 
Y  AND  Ba  SIGNALS. 


FIGURE  6 

COMPUTER  DECONVOLUTION  OF  THE 

INTERFERING  Cu.Y  AND  Ba  SIGNALS 
SHOWN  IN  FI6URE  5. 
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TABLE  3 


RESULTS  FOR  SAMPLE  D 


ELEMENT 

Nt 

REL.  ATOMIC  % 

Y 

3.14+0.15 

1 1 .6+0.3 

Ba 

10.2+0.3 

37.6+0.3 

Cu 

13.8+0.5 

50.8+0.4 

Nt  =  AREAL  DENSITY  (El 7  ATOMS/SQ  CM) 


FIGURE  7 


FITTING  FUNCTION  REPRESENTING 
THE  BACKSCATTERING  SPECTRUM 
OF  A  SINGLE  ELEMENT.  PARAMETERS 
1  AND  2  ARE  WIDTHS  OF  LEFT  AND 
RIGHT  HALF-6AUSSIANS;  PARA¬ 
METERS  3,4,5  AND  6  ARE  HEIGHTS 
AND  LOCATIONS  OF  THE  CENTRAL 
TRAPEZOID. 
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FIGURE  8 

BASIC  BACKSCATTERING 
SPECTROMETRY 

Q  =  NUMBER  OF  HE  IONS  INCIDENT  . 

A*=  COUNTS  IN  BS  PEAK  FOR  ELEMENT  l 
-/2.=  SOLID  ANGLE  SUBT.  BY  DETECTOR . 
GT  =  SCATTERING  XSECTION,  He  ON  {  . 
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